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We construct the field theory which describes the universal properties of the quasi-static isotropic depinning 
transition for interfaces and elastic periodic systems at zero temperature, taking properly into account the non- 
analytic form of the dynamical action. This cures the inability of the 1-loop flow-equations to distinguish 
between statics and quasi-static depinning, and thus to account for the irreversibility of the latter. We prove 
two-loop renormalizability, obtain the 2-loop /3-function and show the generation of "irreversible" anomalous 
terms, originating from the non-analytic nature of the theory, which cause the statics and driven dynamics to 
differ at 2-loop order. We obtain the roughness exponent £ and dynamical exponent z to order e 2 . This allows 
to test several previous conjectures made on the basis of the 1-loop result. First it demonstrates that random- 
field disorder does indeed attract all disorder of shorter range. It also shows that the conjecture C, = e/3 is 
incorrect, and allows to compute the violations, as £ = §(1 + 0.14331e), e = 4 — d. This solves a longstanding 
discrepancy with simulations. For long-range elasticity it yields £ = |(1 + 0.39735e), e = 2 — d (vs. the 
standard prediction £ = 1/3 for d = 1), in reasonable agreement with the most recent simulations. The high 
value of £ ~ 0.5 found in experiments both on the contact line depinning of liquid Helium and on slow crack 
fronts is discussed. 
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I. INTRODUCTION 

A. Overview 

Pinning of coherent structures by quenched disorder, and 
one of its most striking manifestations, the depinning tran- 
sition, are important, ubiquitous and not fully understood 
phenomena 13 . Even a single particle in a quenched random 
potential exhibits a depinning threshold at zero temperature: 
Unbounded motion occurs only when the additional external 
applied force / exceeds a critical force f c . Depinning also 
occurs for systems with many interacting particles, and de- 
pending on the degree of order in the structure, it ranges from 
the so-called plastic depinning 4 to elastic depinning. Here we 
focus on elastic depinning where the particles form a lattice 
or more generally a well ordered structure. The depinning 
transition is then a rather non-trivial collective phenomenon, 
intrinsically out of equilibrium and irreversible: It is well 
known for instance to be a source of hysteresis in magnets 
and superconductors 5 . 

For many experimental systems which exhibit a depinning 
transition a modelization in terms of an elastic object pinned 
by random impurities is a good starting point. The type of 
disorder, which they experience, depends on their symme- 
tries and their local environment. Domain walls in magnets 6 , 
whose study is of importance to information storage technol- 
ogy, behave as elastic interfaces and can experience either 
random-bond disorder (RB), which is short range (SR), or 
random-field disorder (RF), which has long range (LR) spatial 
correlations. Dislocation lines in metals exhibit a depinning 
threshold as the stress in increased 7 . Charge density waves 
(CDW) in solids exhibit a similiar conduction threshold. If the 
applied electric field becomes large enough, the CDW starts 
to slide 8 . Being periodic objects the disorder they feel is also 
periodic 9 . This is also the case for superconductors, where 
vortex lines form, in presence of weak disorder, a quasi or- 



dered periodic Bragg glass phase 10 ' 11 . These systems have 
similarities with (vortex free) continuous XY spins in pres- 
ence of random fields, and generally constitute the random- 
periodic (RP) universality class. 

The contact line of a liquid helium meniscus on a rough 
substrate can be thought of as an interface, but is gov- 
erned by long range elasticity and so are slowly propagating 
cracks 2 ' 1215 . Solid friction is another example of a depinning 
phenomenon. Of course, in each of these systems it must be 
checked separately whether the elastic description holds for 
depinning. It is far from obvious that this is true for all rele- 
vant scales. In any case, in order to be capable to confirm or 
rule out such a description, it is necessary to first obtain pre- 
cise theoretical predictions for the expected behavior in the 
case of elastic depinning, what we aim to achieve here. 

It was proposed some time ago, starting from the study of 
a fully connected mean-field-type model 16 , that the elastic de- 
pinning transition can be viewed in the framework of standard 
critical phenomena. The ordered phase is then the moving 
phase with force / > f c , and the order parameter the veloc- 
ity v, which vanishes as v ~ (/ — f c )@ at the critical point 
/ = f c . The analogy with standard critical phenomena in a 
pure system however has some limits: Additional fluctuation 
exponents were later identified 9 ' 17 , and some non-universality 
was noticed in the fully connected model 1618 . 

It is thus important to develop a renormalization-group de- 
scription of depinning. An important step in that direction 
was performed within the framework of the so-called Func- 
tional Renormalization Group (FRG), to 1-loop order using 
the Wilson scheme 1719-21 . The upper critical dimension was 
identified as d uc = 4, d being the internal dimension of the 
elastic manifold. The peculiarity of the problem is that for 
d < d uc = 4 an infinite set of operators becomes relevant, pa- 
rameterized by a full function A(it), the second cumulant of 
the random pinning force. This problem turns out to be closely 
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related to the statics, i.e. describing the pinned state with mini- 
mal energy in the absence of an applied force / = for which 
the FRG was initially developed 22 (there, the flowing function 
is the second cumulant R(u) of the random potential). Both 
problems are notably difficult due to so-called dimensional 
reduction (DR) which renders the naive T = perturbation 
theory useless 6,23-27 . Indeed to any order in the disorder at 
zero temperature T = 0, any physical observable is found to 
be identical to its (trivial) average in a Gaussian random force 
(Larkin) model. This phenomenon is not restricted to elastic 
manifolds in disorder, but occurs in a broad class of disordered 
systems as e.g. random field spin models and solving it here 
may open the way to a solution in other models as well. The 
FRG at depinning and in the statics seems to provide a way out 
of the DR puzzle: the key feature is that the coarse grained 
disorder correlator becomes non-analytic beyond the Larkin 
scale L c , yielding large-scale results distinct from naive per- 
turbation theory, which assumes an analytic disorder correla- 
tor. Explicit solution of the 1-loop Functional RG equation 
(FRG) for the disorder correlators R(u) and A(u) gives sev- 
eral non-trivial attractive fixed points (FP) 10 ' 22 and critical ex- 
ponents for statics and depinning 10,17 ' 19 ' 21 ' 22 to lowest order in 
e = 4 — d. All these fixed points exhibit a "cusp" singularity, 
which has the form A*(u) — A* (0) <~ \u\ at small The ex- 
istence of the cusp nicely accounts for the existence of a criti- 
cal threshold force 19 , as it is found that f c <~ ^| u _ Q+ A*(u). 

There are however several highly unsatisfactory and puz- 
zling features within the 1-loop treatment, which prompted 
the present and related works. First it was found that the FRG 
flow equation for the statics and depinning are identical to one 
loop (with A(u) = —R"(u)). This implies for instance that, 
within a given universality class (RB,RF and RP), the 1-loop 
RG is a priori unable to distinguish static observables, such as 
the roughness exponent £ at zero applied force / = from 
those at depinning / = f c . This is a rather surprising and 
unphysical result since one knows that depinning is an irre- 
versible out of equilibrium process, quite different from the 
statics. In an attempt to recover the expected physics, and 
to extend conclusions from the 1-loop study to higher orders, 
three conjectures were put forward 17 ' 19-21 : 

1 . At more than 1 -loop order depinning should differ from 
statics. 

2. At depinning the RB universality class should flow to 
the RF universality class: Indeed, since for / — > /+ 
the manifold does not move backward it cannot feel the 
"potential" character of RB disorder. 

3. The roughness exponent of the RF universality class at 
depinning is ( = e/3 to all orders (the Narayan Fisher 
(NF) conjecture 17 ' 21 ), with e = 4 — d for standard man- 
ifold elasticity and e = 2 — d for LR elasticity. 

While conjectures 1 and 2 seem reasonable on physical 
grounds, we emphasize that they were based on qualitative 
arguments: In the absence of any (renormalizable) theory be- 
yond one loop, they appear putative. A 1-loop study including 
the effect of a finite velocity 28 indeed indicated that 2 is cor- 
rect. It strongly relies on a finite velocity, and the behavior in 



the limit v = + was found to be subtle and difficult to fully 
control within that approach. 

The NF conjecture 3 is based on a study of the structure of 
higher orders, but it lacks a controlled field theory argument. 
With the time, it got more and more in disagreement with nu- 
merical simulations and experiments, as we discuss below. In 
addition, if one considers that this value ( = e/3 is expected 
instead for the statics RF class, the NF conjecture seems rather 
unnatural. 

There are also more fundamental reasons to study the FRG 
beyond one loop. In the last fifteen years since 16 ' 22 , no study 
has addressed whether the FRG yields, beyond one loop, a 
renormalizable field theory able to predict universal results. 
There have been 2-loop studies previously but they assumed 
an analytic correlator and thus they only applied below the 
Larkin length 29-31 . Doubts were even raised 32 about the va- 
lidity of the e-expansion beyond order e. 

The aim of the present paper is to develop a more system- 
atic field theoretic description of depinning which extends be- 
yond one loop. A short summary of our study was already 
published 33 together with a companion study on the statics. 
The main and highly non-trivial difficulty is the non-analytic 
nature of the theory (i.e. of the fixed-point action) at T = 0, 
which makes it a priori quite different from conventional crit- 
ical phenomena. It is not even obvious whether this is a legit- 
imate field theory and how to construct it. For the depinning 
transition with / = /+, which is the focus of the present 
paper, we are able to develop a meaningful perturbation the- 
ory in a non-analytic disorder which allows us to show renor- 
malizability at 2-loop order. Even the way renormalizability 
works here is slightly different from the conventional one. To 
handle the non-analyticity in the static problem is even more 
challenging, and we propose a solution of the problem to 2- 
loop 33 34 and 3-loop 35 order as well as at large- N 36 . 

In this paper we focus on the so-called "isotropic depin- 
ning" universality class. This means that the starting model 
has sufficient rotational invariance, as discussed below, which 
guarantees that additional Kardar-Parisi-Zhang terms are ab- 
sent. A general discussion of the various universality classes 
can be found in 37 38 and an application of our non-analytic 
field theory (NAFT) methods to the case of "anisotropic de- 
pinning" will be presented in 39 . 

Before we summarize the novel results of the present pa- 
per, let us recall some important features about the model, the 
scaling and statistical fluctuations at the depinning threshold. 



B. Model, scaling and fluctuations 

Elastic objects can be parameterized by a TV-component 
height or displacement field u x , where x denotes the d- 
dimensional internal coordinate of the elastic object (we will 
use u q to denote Fourier components). An interface in the 3D 
random field Ising model has d = 2, N = 1, a vortex lat- 
tice d = 3, N — 2, a contact-line d = 1 and N = 1. In 
this paper we restrict our study to N = 1. In the presence of 
a random potential the equilibrium problem is defined by the 
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Hamiltonian: 

H = J^-u q u- q + jv{u x ,x) (1.1) 

with c(q) = cq 2 for standard short-range elasticity, c(q) = 
c\q | for long-range elasticity and we denote J q = J and 
J = J d d x. Long-range elasticity appears e.g. for the contact 
line by integrating out the bulk-degrees of freedom 40 . For pe- 
riodic systems the integration is over the first Brillouin zone. 
More generally a short scale UV cutoff is implied at q <~ A, 
and the system size is denoted by L. As will become clear 
later, the random potential can without loss of generality be 
chosen Gaussian with second cumulant 



V(u, x)V(u', x>) = R(u - u')S d (x - x) . (1.2) 

Periodic systems are described by a periodic function R(u), 
random bond disorder by a short range function and random 
field disorder of amplitude a by R(u) <~ — a\u\ at large u. 

We study the over-damped dynamics of the manifold in this 
random potential, described (in the case of SR-elasticity) by 
the equation of motion 

ridtu x t = c\7 2 x u xt + F(x, u xt ) + f (1.3) 

with friction r\. In presence of an applied force / the center 
of mass velocity is v = L~ d J dtu x t- The pinning force is 
F(u, x) — — d u V(u, x) and thus the second cumulant of the 
force is 



F(x, u)F(x\ u') = A(u - u')5 d (x - x') , (1.4) 

such that A(u) = —R"(u) in the bare model. As we will 
see below it does not remain so in the driven dynamics. The 
"isotropic depinning" class contains more general equations 
of motion than (1.3). For instance some cellular automaton 
models are believed to be in this class 41 . They must obey 
rotational invariance, as discussed in Ref. 37 ~ 39 , which prevents 
the additional KPZ term A(V ' x u xt ) 2 to be generated at / = 
/+. There is always a KPZ term generated at v > from 
the broken symmetry x — ► —x, but A can vanish or not as 
v — > + , depending on whether rotational invariance is broken 
or not. Here this symmetry is implied by the statistical tilt 
symmetry (STS) 42 ' 43 u xt — > u xt + g x . It also holds in the 
statics and accounts for the non-renormalization of the elastic 
coefficient, here set to c = 1. 

A quantity measured in numerical simulations and exper- 
iments is the roughness exponent at the depinning threshold 

f = fc 

C L (x - x 1 ) = \u(x) - u{x')\ 2 ~ \x - x'\ 2i , (1.5) 

which can be compared to the static one (cq- Other expo- 
nents have been introduced 1617 ' 19-21 . The velocity near the 
depinning threshold behaves as«~ (/ — fc) 13 ', the dynamical 
response scales with the dynamical exponent t ~ x z and the 
local velocity correlation length £ diverges at threshold with 
£ ~ (/ — Sc) v - There have also been some studies below 



threshold 9 44 . The following exponent relations were found to 
hold 19 : 

P = v(z- C) (1-6) 

v - 2^ (L7) 

the latter using STS. There are various ways to measure the 
roughness exponent. In some simulations 45 ^ 17 it has been ex- 
tracted from the critical configuration, i.e. as / is increased to 
f c in a given sample it is obtained from the last blocking con- 
figuration. It can also be defined as the limit v — > + of the 
roughness in the moving state, which we will refer to as the 
"quasi-static" depinning limit to distinguish it from the previ- 
ous one. This is the situation studied in this paper. Although 
it is widely believed that both are the same, the depinning the- 
ory has enough peculiarities that one should be careful. In 
particular, beyond scaling arguments and simulations, there is 
presently no rigorous method capable to connect the behavior 
below and above threshold. 

Another peculiarity was noted in 17 . It was found that the 
finite-size fluctuations of the critical force can scale with a 
different exponent: 

f c (L) - U - L~ lf " FS (1-8) 
and it was questioned whether i^ps = v. The bound 

^ FS >2/(rf + C) (1-9) 

follows from a general argument of Ref. 48 . For charge density 
waves where £ = one sees that v = 1/2 and thus v and 
t'FS must be different for d < 4. For interfaces it was noted 17 
that v = z/ps is possible provided £ > e/3. If one assumes 
v = fps, the NF-conjecture £ = e/3 is then equivalent to 
saturating the bound (1.9). We will address the question of 
whether v = i/ps below. 

Finally note that at / = f c the condition of equilibrium of 
a piece of interface expresses that the elastic force, which acts 
only on the perimeter, balances the excess force on the bulk, 
yielding the scaling: 

L d - l u{a,L)~{f c {L)- f c )L d , (1. 10) 

where u(a, L) <~ \/CiJa) is the relative displacement (1.5) 
between two neighbors averaged over the perimeter. This 
shows that 

u(a,L) - L 1 ^^ (I.H) 

thus for CDW the displacements between two neighbors 
grows unboundedly 49 with L for d < 2. For interfaces (non- 
periodic disorder), if one assumes v — ^ FS one obtains that 
the displacements between two neighbors grows with L only 
when £ > 1. 

C. Summary of results 

Let us now discuss the main results of our study. 
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First we show that, at depinning, 1- and 2-loop diagrams 
can be computed using a non-analytic action in an unambigu- 
ous and well defined way, allowing to escape dimensional re- 
duction. The mechanism is non-trivial and works because the 
manifold only moves forward in the steady state which allows 
to remove all ambiguities. We show that the limit v — > + 
can be taken safely without additional unexpected singulari- 
ties arising in this limit. 

Next we identify the divergences in the 2-loop diagrams us- 
ing dimensional regularization in d = 4 — e. We identify the 
1 -loop and 2-loop counter-terms and perform the renormaliza- 
tion program. We find that the 1/e divergences cancel nicely 
in the /3-function for the disorder correlator and in the dynam- 
ical exponent. The theory is finite to 2-loop order and yields 
universal results. 

The obtained FRG flow equation for the disorder (the j3- 
function) contains new "anomalous" terms, absent in an an- 
alytic theory (e.g. in the flow obtained in Ref. 29 ' 30 ). These 
terms are different in the static theory (obtained in 33 ) and at 
depinning, showing that indeed static and depinning differ to 
two loops. Thus the minimal consistent theory for depin- 
ning requires two loops. 

Next we study the fixed point solutions of our 2-loop FRG 
equations at depinning. For non-periodic disorder (e.g. inter- 
faces) with correlator of range shorter or equal to random- 
field, we find that there is a single universality class, the 
random-field class. Thus random-bond disorder does flow to 
random field. Specifically we find that the flow of J A is cor- 
rected to two loops and thus J A cannot remain at its random- 
bond value, which is zero. This is explained in more detail in 
section IV. The problem does not remain potential and irre- 
versibility is manifest. For short range elasticity, we find the 
roughness-exponent at depinning: 

C = |(l + 0.143313 e) (1.12) 

with e = 4 — d, and for long range elasticity: 

C= |(l + 0.39735 e) (1.13) 

with e = 2 — d. Thus the NF-conjecture 17 ' 21 that £ = | is 
incorrect. We also compute the dynamical exponent z and 
obtain (3 and v by the scaling relations (1.6) and (1.7). We also 
find that vp?, = v holds to two loops. 

For periodic disorder, relevant for charge density waves, 
we find a fixed point which leads to a universal logarithmic 
growth of displacements. This fixed point is however unsta- 
ble, as an additional Larkin random force is generated. The 
true correlations are the sum of this logarithmic growth and of 
a power law growth so that the true ( = (4 — d)/2. This is 
similar to 50 . Then we find 

v = \ (1-14) 
^fs = 3 . (1-15) 



D. Numerical simulations and experiments 

Over many years, numerous simulations near 
depinning 19 ' 20 ' 5 1-54 accumulated evidence that £ 7^ e/3. 
In d = 1 in particular often an exponent ( > 1 was observed. 
Our results show that ( > e/3 and thus resolve this long 
standing discrepancy between numerical simulations and the 
renormalization group. They are summarized in Tables IV.2 
and IV.3 in Section IV, where we compare them to numerical 
simulations. Of course it is not possible to give strict error 
bars from the FRG calculation without further knowledge of 
higher orders, but one can still give rough estimates, based on 
different Pade-approximants. 

Let us in the following discuss recent numerical results. 
Following shortly our paper 33 , Rosso and Krauth obtained a 
set of precision numerical results using a powerful algorithm 
to determine the critical configuration at depinning (the last 
blocking configuration) up to large sizes 45 "* 7 . They obtained 
results in d = 1 which, despite being far from d = 4, compare 
well with our results. For short range elasticity they find 

C = 1.17 (1.16) 

close to our 2-loop result (1.12). Note that displacement cor- 
relations scaling as 

ujn- q ~ <r (d+2 ° (1-17) 
with Q > 1 are perfectly legitimate. It simply means that 

C L (x) - 2 f(l - cos(qx))q-( d+2 V ~ L 2( <-^x 2 . (1.18) 

Jq 

The size dependent factor comes from the infrared divergence 
of the integral. Thus in a simulation neighboring monomers 
will be spread further and further apart, which is fine if their 
attraction is purely quadratic. Of course in a realistic physical 
situation their bond will eventually break, but as a model it 
is mathematically well defined. For the anisotropic depinning 
universality class, not studied here, they found ( = 0.63 as 
many other authors using cellular automaton models 55-57 . 

For isotropic depinning with long range elasticity they ob- 
tained: 

C = 0.390 ± 0.002, (1.19) 

which lies roughly at midpoint of the 1-loop and 2-loop pre- 
diction setting e = 1 in (1.13). So do their most recent 
estimates 58 forSR disorder. In d = 2 this is ( = 0.753+0.002 
and for d = 3 they obtain 0.35 < ( < 0.4. These results (1.16) 
and (1.19) are close to estimates from the 2-loop expansion 
and clearly rule out the NF-conjecture. 

Another recent work 59 studies an interface in the random 
field Ising model in high dimension. The authors confirm that 
d = 4 is the upper critical dimension. They further extract the 
velocity exponent (3 and compare their results with our 2-loop 
FRG prediction for (3: 



which holds presumably to all orders. 
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(1.20) 
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FIG. 1.1: Figure from Ref. which compares new numerical values 
(black circles) and a previous one (white square) obtained for the 
exponent (3 with our prediction from the FRG. 

The results are shown in Fig. 1. 1 . One can see a clear curvature 
downwards and that the straight line giving the 1-loop result 
is well above the obtained results (the 1-loop approximation 
would predict (3 = 0.78 in d = 2). 

Thus, although there is still some spread and uncertainty 
in the results, it seems that there is now a trend towards a 
convergence between theory and numerical simulations. 

The situation concerning experiments is presently unclear. 
Let us first outline the generic findings before analyzing 
the details. The measured exponents corresponding to LR- 
elasticity and d = 1 seem to be consistently in the range 
C ~ 0.5 — 0.55. This is slightly above our 2-loop result (1.13) 
but not fully incompatible with it. Our calculation holds for 
quasi-static depinning, i.e. v > — > + , and most exper- 
iments are also performed from the moving side, hopefully 
reaching the same quasi-static limit v — > + . On the other 
hand if one believes that the numerical result 1.19 (also com- 
patible with our calculation, from below) obtained for / = f~ 
also holds for quasi-static depinning (a rather natural, but as 
yet unproved assumption) then one must conclude that the 
elastic models, in their simplest form at least, may not faith- 
fully represent the experimental situation. Care must how- 
ever be exercised before any such conclusion is reached. One 
could argue that disorder A(u) <~ u~ a of range longer than 
RF (a < 1) could produce higher exponents £ = e/(2 + a) 
(see end of Section IV A) but that does not seem to apply to 
those experiments where disorder is well controlled. Also, 
since the exponent ( = 0.5 is the Larkin DR-exponent, which 
should hold below the Larkin length L c one must make sure 
that L c is well identified and that one is not simply observing 
a slow crossover to the asymptotic regime. In some of these 
experiments L c has been identified to be rather small. 

Let us now examine the situation in more details. 

One much studied experimental system is the contact line 
of a fluid 12 60 . It advances on a rough substrate and is pushed 



by adding fluid to the reservoir. The elasticity of the line is 
short range at short scale but at larger scales it is mediated 
by the elasticity of the two dimensional meniscus and thus 
it becomes long range and should be compared with (1.13), 
(1.19). Disorder is random-field, but one should distinguish 
between microscopic disorder, which is poorly characterized, 
and macroscopic one which is well controlled. The situation 
has been studied for a helium meniscus on a macroscopically 
disordered substrate where £ = 0.55 was found 60 . Although 
there are good indications that these experiments probe quasi- 
static depinning (the contact line jumps from a reproducible 
pinned configuration to the next one) the precise nature of the 
dynamics remains open. Indeed it was found that propagation 
of perturbations along the line can be as fast as avalanches, 
showing inertial regime for helium 68 . Experiments were re- 
peated for viscous liquids 69 yielding ( = 0.51 ± 0.03. There 
it was checked that the system is over-damped and near depin- 
ning. In both cases there is also evidence of thermal activation 
effects 13 characteristic of depinning (not creep). It was argued 
that these may be a signature that a more complicated dynam- 
ics (e.g. plastic) takes place at the very short scales and pro- 
duces an effective dynamics at larger scales with complicated 
non-linear (e.g. exponential) velocity and temperature depen- 
dent damping. Very similar effects have also been shown to 
occur in solid friction 70 were the activation volume was also 
found to correspond to microscopic scales. 

Another class of much studied experimental systems are 
crack fronts in heterogeneous media 71 . These are charac- 
terized by two displacement fields, one out-of-plane compo- 
nent h and an in-plane one /. Cracks can either be studied 
stopped or slowly advancing. At the simplest level the in- 
plane displacement / is expected to be described as an elas- 
tic line d = 1, N = 1 with LR elasticity c\q\, at quasi 
static depinning 72 . In experiments 15,73 the observed rough- 
ness is again Q 0.55. Since the crack propagates in an 
elastic medium, elastic waves which can in principle affect 
the roughness as the crack front advances producing a more 
complicated dynamics than Eq. (1.3). Some proposals have 
been put forward on mechanisms to produce higher roughness 
exponents 74 They rely however on a finite velocity and it is 
unclear whether they can modify roughness in the quasi-static 
limit. Even if instantaneous velocities during avalanches be- 
come large enough, a detailed description on how these could 
change the line configurations remains to be understood. Then 
of course a major issue is whether the experiment, and in 
which sense, is in the quasi-static limit. There again micro- 
scopic dynamics could be more complex as at small scales 
the material may be damaged and the notion of a single front 
may not apply. Finally, since there are two components to 
displacement one should also be careful to understand inter- 
actions between them near depinning 61 . 

Another interesting experimental system is a domain wall 
in a very thin magnetic film 62 which experiences RB disorder. 
Up to now however only the thermally activated motion has 
been studied, which gives a quite remarkable confirmation of 
the creep law 62 with RB exponents. It would be interesting 
to study depinning there and to check whether it also belongs 
to the isotropic universality class. In that case, the crossover 
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from RB to RF resulting in overhangs beyond some scale at 
zero temperature (( > 1) as well as the non-trivial thermal 
rounding of depinning could be studied. 

II. MODEL AND PERTURBATION THEORY 

In this Section we discuss some general features of the field 
theory of elastic manifolds in a random potential, both for the 
statics and for the dynamics, driven or at zero applied force. 
Some issues are indeed common to these three cases. At the 
end we specialize to depinning. 

A. Static and Dynamical action and naive power 
counting 

The static, equilibrium problem, can be studied using replicas. 
The replicated Hamiltonian corresponding to (1.1) is: 

Jx a 

where, for now, we consider SR elasticity, a runs from 1 to n. 
We have added a small mass to provide an infrared cutoff, and 
we are interested in the large scale limit m — > 0. The limit of 
zero number of replicas n = is implicit everywhere. Terms 
with sums over three replicas or more corresponding to third 
or higher cumulants of disorder are generated in the perturba- 
tion expansion. These should in principle be included, but as 
we will see below higher disorder cumulants are not relevant 
for the T = depinning studied below. 

The dynamics, corresponding to the equation of motion 
(1.3) is studied using the dynamical action averaged over dis- 
order: 

S[u,u] = / iu xt {r\d t - <9 2 + m 2 )u xt - rjT iu xt iu xt 

J xt J xt 

-- [ iuxtiuxt' M u xt - u x f) - / iiixtfxt ■ 

* Jxtt' Jxt 



(11.2) 



It generates disorder averaged correlations, e.g. (A[u xt }} — 
(A[u xt ])s with (A) s = J V[u]V [u]Ac- s and (l) s = 1, and 
response functions 5(A [u])/5f x t — (iu x tA[u])s- The uni- 
form driving force f xt — f > (beyond threshold at T = 0) 
may produce a velocity v = dt{u x t) > 0, a situation which 
we study by going to the comoving frame (where (u xt ) = 0) 
shifting u xt — ► u xt + vt, resulting in f —* f — r/v. This is 
implied below. In general, for any value of /, we study the 
steady state, which at finite temperature T > is expected 
to be unique and time translational invariant (TTI) (all aver- 
ages depend only on time differences). In the zero tempera- 
ture limit, one needs a priori to distinguish the T = TTI 
theory as Hin^^oo linir^o (e.g. the ground state in the static) 
and the T = + theory as liniT^o liniL^oo- 

It is important to note that there are close connections, via 
the fluctuation dissipation relations, between the dynamical 



formalism and the statics. Indeed, at equilibrium (for / = 
and when time translation invariance is established) any equal 
time correlation function computed with (II.2) is formally 
identical (e.g. to all orders in perturbation theory) to the corre- 
sponding quantity computed in the equilibrium theory (which 
is a single replica average). Similarly, the persistent parts, i.e. 
those oc S(oj), of dynamical correlations involving p mutually 
very separated times, are formally identical to the correspond- 
ing averages in the replica theory involving p replicas. The 
perturbative equilibrium calculations in the statics can thus be 
indifferently performed either with replicas or with (II.2). It 
is possible to generate all dynamical graphs from static ones, 
a connection which, as will be further explained below, also 
carries to some extent to the case / > at T = 0. 

We first study "naive" perturbation theory and power count- 
ing. The quadratic part So of the action (II. 2) yields the free 
response and correlation functions, used for perturbation the- 
ory in the disorder. They read 



\ D 0(t-t') _I*=£l ( „2 +m 2 } 

{iUq t t>u- q>t )o = R q .t-t' = — -e " {q + ; 



(0.3) 



respectively, with the FDT relation TR q T — —d T C q . T (r > 
0). Perturbation theory in A(u) yields a disorder interaction 
vertex and at each (unsplitted) vertex there is one conservation 
rule for momentum and two for frequency. It is thus conve- 
nient to use splitted vertices, as represented in Fig II. 1, where 
the rules for the perturbation theory of the statics using replica 
are also given. For the dynamics one can also focus on T = 
where graphs are made only with response functions and con- 
sider temperature as an interaction vertex. The 1-loop and 2- 
loop diagrams which correct the disorder at T = are shown 
in Fig. D.l (unsplitted vertices). There are three types of 2- 
loop graphs A, B, C. The graphs E and F lead to corrections 
proportional to temperature. 

At T — the model exhibits the property of dimensional 
reduction 6,23-27 (DR) both in the statics and dynamics. Its 
"naive" perturbation theory, obtained by taking for the dis- 
order correlator A(u) an analytic function of u (or R(u) 
for the statics), has a triviality property. As is easy to show 
using the above diagrammatic rules the perturbative expan- 
sion of any correlation function (]J i u Xiti )s in the derivatives 
A(")(0) yields to all orders the same result as that obtained 
from the Gaussian theory setting A(m) = A(0) (the so called 
Larkin random force model). The same result holds for the 
statics, for any correlation (]J t u%\)s- At T = these correla- 
tions are independent of the replica indices their dynamical 
equivalent being independent of the times U. The two point 
function thus reads to all orders: 



(u q ,t'U- q ,t) BR = 



A(0) 



(g 2 +m 2 ) 



2\2 



(0.4) 



This dimensional reduction results in a roughness exponent 
£ = (4 — d)/2 which is well known to be incorrect. One 
physical reason, in the statics, is that this amounts to solv- 
ing the zero force equation which, whenever more than one 
solution exists, is not identical to finding the lowest energy 
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FIG. II. 1: (i) diagrammatic rules for the statics: replica propaga- 
tor (u a Ub)o = TSab/q 2 , unsplitted vertex, equivalent splitted ver- 
tex — 5~Jaf> TfiTiRiua — ?*b) and (ii) dynamics: response propagator 
(ww)o = R q ,t-t', unsplitted vertex, splitted vertex u xt u xt i A(u xt — 
u xt i ) and temperature vertex. Arrows are along increasing time. An 
arbitrary number of lines can enter these functional vertices, (iii) 
unsplitted diagrams to one loop D, one loop with inserted 1-loop 
counter-term G and 2-loop A,B,C,E,F. 

configuration. Curing this problem, within a field theory, is 
highly non-trivial. One way to do that, as discussed later will 
be to consider a non-analytic A(u). 

It is important to note that despite the DR, dynamical aver- 
ages involving response fields remain non-trivial, even at zero 
temperature. Perturbation theory at finite temperature also re- 
mains non-trivial. It is thus still useful to do power counting 
with an analytic A(u), the modifications for a non-analytic 
A(m) being discussed in the following section. 

Power counting at the Gaussian fixed point yields t <~ x 2 
and uu <~ x~ d . At T = nothing else fixes the dimensions 
of u, since u — > Am, u — ► \~ 1 u leaves the T = action 
invariant. Denoting u <~ x^, ( is for now undetermined. The 
disorder term then scales as x 4 ~ d+2 ^. It becomes relevant for 
d < 4 provided £ < (4 — d)/2 which is physically expected 
(for instance in the random periodic case, ( — is the only 
possible choice, and for other cases £ = 0(e)). With this 
power counting the temperature term scales as x~ e with 9 = 
d—2 + 2( and is thus formally irrelevant near four dimension. 
In the end £ will be fixed by the disorder distribution at the 
fixed point 22 . 

A more detailed study of divergences in the vertex func- 
tions allows to identify all counter-terms needed to render the 
theory finite. We denote by 

Fu..u;u..u(qi,Ui,qi,Ui) = (H.5) 

n ^ — n m~ t[u ' fi] i«=*=o 

the irreducible vertex functions (IVF) with E u external fields 
u (at momenta q t ,uii,i = 1, ..E u ) and E u external fields (at 
momenta qi,u>i, i = 1, ..En). Being the derivative of the ef- 
fective action functional T[u, u] they are the important objects 
since all averages of products of u and u fields are expressed 
as tree diagrams of the IVF. Finiteness of the IVF thus im- 
ply finiteness of all such averages. The present theory has 




O 




FIG. II.2: Construction of diagrams starting from an unsplitted static 
diagram via two splitted static diagrams (2-replica component) to the 
corresponding dynamical diagrams as explained in the text. 

the property of covariance under the well known statistical tilt 
symmetry STS u xt — > u xt + g x , which yields that the two 
point vertex Tu U (uj = 0) remains uncorrected to all orders. 
This allows to fix the elastic constant c = 1 and shows that 
the mass term is uncorrected and can thus safely be used as an 
IR cutoff. It also implies that all higher IVFs vanish when any 
of the oji is set to zero. The DR result is a perturbative triv- 
iality statement about T^^uiQi, &i) at T = 0, all other cases 
remain non-trivial. In a sense we will now expand around di- 
mensional reduction. Similar replica IVFs can be defined for 
the statics. 

Perturbation expansion of a given IVF to any given order in 
the disorder can be represented by a set of one particle irre- 
ducible (1PI) graphs. As mentionned above there is a simple 
rule to generate the dynamical graphs from the static ones. 
The static propagator being diagonal in replicas, each static 
graph occurring in a p replica IVF contains p connected com- 
ponents. At T = the rule is then to attach one response 
field to each connected component of the static diagram, each 
replica graph then generating one or more dynamical graphs. 
The place where the response field is attached is the root of 
the diagram. The direction of the remaining response func- 
tions is then fixed unambiguously, always pointing towards 
the root. This procedure to deduce the dynamical diagrams 
from the static ones is unique and exhaustive and is illustrated 
in Fig. II. 2. A generalization exists at T > but is not needed 
here. 

Any graph corresponding to a given dynamical IVF con- 
tains p connected components (in the splitted diagrammatics) 
with 1 < p < Eu (p = Eu at T = 0), each one leading to a 
conservation rule between external frequencies, and thus one 
can write symbolically: 

^u..u-,u..u{qii Qi: ^i) 

= *(E^-E?)rKE--E-)r. ( n - 6 > 

i=l 

Let us compute the superficial degree of UV divergence a of 
such a graph with va disorder vertices and vt temperature 
factors contributing to T ~ A". Using momentum and fre- 
quency conservation laws at each vertex, and since there are 
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only response functions E u + I = 2(va + vt) we obtain: 

a = d + 2p-dE u + (d- 4)v A + {d- 2)v T . (11.7) 

At T — (vt = 0, p = E u ) at the critical dimension d = 4 
the only superficially UV divergent IVF are those with one ex- 
ternal u (quadratic divergence) or two external u (logarithmic 
divergence LD). The STS further restricts the possible diver- 
gent diagrams. One sees that only three types of counter-term 
are needed a priori. One counter-terms is needed for the A 2 
divergence of T u (q = 0, u = 0) (excess force / — i]v in driven 
dynamics). This is analogous to the mass in the </> 4 theory, i.e. 
the distance to criticality. If we are exactly at the depinning 
critical point (/ = f c ) we need not worry about this diver- 
gence. Another counter-term is associated with the LD in r] 
and the last one with the LD in the second cumulant of disor- 
der A(u), i.e. a full function, which makes it different from the 
conventional FT for critical phenomena (e.g. </> 4 ). One notes 
that higher cumulants are formally irrelevant, as they involve 
Eu > 2. 

One sees from (II. 7) that each insertion of a temperature 
vertex yields an additional quadratic divergence in d = 4, 
more generally a factor TA d ~ 2 . Thus to obtain a theory where 
observables are finite as A — > oo one must start from a model 
where the initial temperature scales with the UV cutoff as 



T = Tm 



d-2 



(H.8) 



This is similar to the cf) 4 theory where it is known that a <j) 6 
term can be present and yields a finite UV limit (i.e. does not 
spoil renormalizability) only if it has the form g e (j) 6 /A d ~ 2 . It 
then produces only a finite shift to 174 without changing uni- 
versal properties 75 . Here each T factor will thus come with a 
A 2 ~ d factor which compensates the UV divergence. Comput- 
ing the resulting shift in A(u) to order A 2 by resumming the 
diagrams E and F of Fig. II. 1 and all similar diagrams to any 
number of loops has not been attempted here. 

For convenience we have inserted factors of m in the defini- 
tion of the rescaled temperature, using the freedom to rescale 
u by m _< » and u by mf. The disorder term then reads as 
in (II. 2) with A(u) replaced by A (u) = to £ ~ 2? A(wm£) in 
terms of a dimensionless rescaled function A. 



B. Non-analytic field theory and depinning in the 
quasi-static limit 

From now on we study the zero temperature limit T = 0. 
To escape the DR triviality phenomenon, and since the fixed 
points found in 1-loop studies exhibit a cusp at u = 0, we 
must consider perturbation theory in a non-analytic disorder 
correlator. In this section we show how to develop perturba- 
tion theory and diagrammatics in a non-analytic theory and 
what are the non-trivial issues which arise. 

For now the considerations apply for zero or finite applied 
force. In usual diagrammatics, extracting a leg from a vertex 
corresponds to a derivation. Here this can be done as usual 
with no ambiguity, provided the corresponding vertex is eval- 
uated at a generic u (e.g. the graphs in Fig II. 2). If the vertex is 



evaluated at u = (here and in the following we call them sat- 
urated vertices) one must go back to a careful application of 
Wick's rules. Any graph containing such a vertex and which 
vanishes in the analytic theory is called anomalous. Let us 
write the series expansion in powers of |u|: 



A{u) = A(0) + A'(0+)M + -A"(0+)u 2 



+ 



(H.9) 



Wick's rules can then be applied but usually end up in evalu- 
ating non-trivial averages of e.g. sign or delta functions. 

Let us consider as an example the following 2-loop 1PI di- 
agram (noted C! in what follows) which is a correction to the 
effective action of the form: 




j-A (Wr.r ^r.cr) 



(11.10) 



x / R ri — r2i t 1 R ri — r2i t 2 R r — ri j 3 R r — r2 j 4 

Jti>0, ri 

xA'(ll ri!T -t 3 — U ri)T _t 4 _t 1 )A'('U T . 2ir _t 4 — u r 2 ,T-t 2 -t 3 ) ■ 

Here four Wick contractions have been performed, as in any 
of the other thirty 2-loop diagrams of the form A (studied in 
the next Section). In an analytic theory performing the lo- 
cal time expansion this would result in a 2-loop correction 
to A(u) proportional to A"(u) but with a zero coefficient 
since the A' functions are evaluated at zero argument. In 
the non-analytic theory, inserting the expansion (II. 9) yields 
(upon some change of variables): 

ei = A'(0+) 2 A'» / R rutl R ri , t2 R r3 

Jti>0,ri 

F riM = (sgn(X)sgn(F)) (11.11) 
A = u rii -t 3 — u ri ^ ti -t 1 
Y = u ,-u - uo,-t 3 -t 2 

terms of higher order in (II. 9) do not contribute since we are at 
T = and we have exhausted the number of u to contract (i.e 
those terms would yield higher orders in T). The remaining 
average in (II. 1 1) is evaluated with respect to a Gaussian mea- 
sure, and can thus be performed. It can be defined by using the 
T > 0, v > Gaussian measure (u xt — > vt + u xt ) and taking 
the limit T — > 0, v — > 0. The result is a continuous function 
ofv 2 /T and its value depends on how the limit is taken. 

In the static theory one should take T — > at v = 0. This 
yields 



(sgn(A)sgn(F)) = — asin(cr) 
(XY) 



(0.12) 



a = 



i.e. the result for centered Gaussian variables. Expressing the 
averages in (11.12) using correlation functions C q _ t yields a 
complicated T = expression for e\. This expression will be 
discussed in a companion paper on the statics 34 . A list of all 
anomalous diagrams is presented in Appendix K. 
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The opposite limit v — > at T = yields much simpler 
expressions: 

(sgn(X)sgn(F)) -> sgn(t 4 + *i - t 3 )sgn(t 3 + i 2 - *4) • 

More generally this procedure corresponds to the substitution: 
A(")(u r ,t — u r ,t> ) A^ n \v(t — t')) in any ambiguous vertex 
evaluated at u = 0. That this is the correct definition of the 
theory of the quasi-static depinning as the limit v = + is 
particularly clear here since it is well known (the no passing 
property 9 ' 44 ) that the u r ^ are increasing functions of time in 
the steady state. Of course it remains to be shown that the 
procedure actually works and does not produce singular terms 
such as 6(vt). It also remains to be shown that it yields a 
renormalizable continuum theory where all divergences can 
be removed by the appropriate counter- terms. This is far from 
trivial and will be achieved below. 

Let us comment again on the connections between dynam- 
ics and statics. Consider a T = dynamical diagram with p 
connected components evaluated at zero external frequencies. 
All response functions can be integrated over the times from 
the leaves towards the root on each connected component. Us- 
ing the FDT relation this replaces response by correlations and 
thus exactly reproduces a p replica static diagram except that 
it is differentiated once with respect to each replica field (the 
sums over all possible positions of the response field repro- 
duces the derivation chain rule). One simple way to establish 
this rule is to consider the formal limit rj — > + (equivalently 
expansion of R QjW in powers of frequency), i.e. formally re- 
place R q ,t,t> — > Su'/q 2 (keeping track of causality). This 
reproduces exactly the zero frequencies dynamical diagrams 
and treats "replicas" as "times". 

Thus the p-th derivative of a p replica static diagram gives 
a set of dynamical diagrams with p connected components. 
Forp = 2 this ensures e.g. that the relation A(u) = —R"(u) 
remains uncorrected to all orders. The flaw in this argument 
comes from the anomalous diagrams (both in statics and dy- 
namics). In the analytic theory the dynamical diagrams with 
response fields on a saturated vertex vanish or cancel in pairs. 
This just expresses that taking a derivative of a static saturated 
vertex gives zero and the rule still works. But in the non- 
analytic theory the anomalous diagrams do not vanish and 
contain an additional time dependence. The above integra- 
tion of response functions from the leaves to the root cannot 
be performed for these anomalous diagrams. As a result they 
can give non-trivial contributions both in statics and dynamics 
which violate relations such as A(u) = —R"(u), thus allow- 
ing to distinguish statics from depinning. 

To conclude this Section: The perturbative calculation of 
the effective action and of the IVF vertices can also be per- 
formed in a non-analytic theory. It can be expressed as sums 
of the same diagrams one writes in the analytic theory, with 
the same graphical rules to draw and generate the diagrams 
starting from the statics. However the way to compute these 
diagrams and their values is different from the analytic theory. 
The time ordering of vertices comes in a non-trivial way and 
produces results which can be different at depinning / = /+ 
(v = + ) and in the statics / = 0, as illustrated on the diagram 
ei above. Thus we see the principle mechanism by which 



the statics and the depinning can yield different field theories, 
which is a novel result. It remains to perform the actual calcu- 
lation of these non-analytic diagrams, which is performed in 
the following sections. 

III. RENORMALIZATION PROGRAM 

In this section we will compute the effective action to 2-loop 
order at T = for depinning. From the above analysis we 
know that we only need to compute the 1- and 2-loop correc- 
tions to A(u) and r\. 

A. Corrections to disorder 

We start by the corrections to the disorder, first at 1-loop and 
then at 2-loop order. 

1. One loop 




FIG. III. 1 : 1-loop dynamical diagrams correcting A 

At leading order, there are four diagrams, depicted on fig- 
ure III. 1 . Since diagram (d) is proportional to A'(u)A'(0), 
it is an odd function of u, and thus does not contribute to 
the renormalization of A. However its repeated counter-term 
will appear at 2-loop order. Diagram (a) is proportional to 
-A(u)A"(u), diagram (b) to -A'(u) 2 and diagram (c) to 
A"(u)A(0). All come with a combinatorial factor of 1/2! 
from Taylor-expanding the exponential function, 1/2 from the 
action and 4 from combinatorics. Together, they add up to the 
1-loop correction to disorder 

^A(n) = A [- A '(uf (A(«) - A(0))A»] h 

J i : = / ( 2 t 2,2 cnu) 

J q (q 2 + m 2 ) 2 
with h = j q c-" 2 r(2 -!) = (47r)- d / 2 r (2-|). 

2. Two loops 

First, we have to find all diagrams correcting disorder at sec- 
ond order. At T = they can be grouped in 3 classes A, B and 
C for the 3 possible diagrams for unsplitted vertices. Class C 
does not contribute as is shown in appendix C. We begin our 
analysis with class A. 




FIG. III. 2: The 3 possible classes at second order correcting disorder 
at T — 0. Only classes A and B will contribute. 




FIG. III.3: Static graphs at 2-loop order in the form of a hat (class A 
in figure III.2) contributing to two replica terms. Adding a response- 
field to each connected component leads to the dynamic diagrams of 
figure ni.4. 

We now need to write all possible diagrams with splitted 
vertices of type A. A systematic procedure is to start from 
all possible static diagrams given in Fig. III. 3. This relies 
on the fact that dynamics and statics are related - recall that 
in general a dynamic formulation can be used to obtain the 
renormalization of the statics. As mentionned in the previous 
Section, to go from the statics to the dynamics, one attaches 
one response field to a root on each connected component of 
the diagrams a to f in figure III. 3 and orient each component 
towards the root. The result is presented on figure III.4. The 
next step is to eliminate all diagrams which yield odd func- 
tions of u and thus do not contribute to the renormalized dis- 
order. The list is the following: 

ai — 8.4 — c 3 = di = d 3 = ds = d7 

= e 2 = e 3 = fi = f 3 = h = {5 = • (HI.2) 

Further simplifications come from diagrams, which mutually 
cancel. Again this uses that A'(u) is an odd function. This 
gives: 

c 2 + c 5 = d 2 + d 4 = d 6 + d 8 = . (III.3) 

In addition 

c 4 = , (III.4) 

since J.., R x tR x t' A' (t — t') = 0. This is explained in more 
details in Appendix K where the list of all anomalous (non- 
odd) graphs is given together with their expressions in the 
non-analytic field theory. 

Thus, the only non-zero graphs which we have to calcu- 
late are a 2 , a 3 , bi,...,b 6 , Ci, ei and f 2 . These calcula- 
tions are rather cumbersome, due to the appearance of theta- 
functions of sums or differences of times as a result of the 
non-analyticity of the theory. The correction to disorder is 

i 2 A( u ) = l|3(2 3 )^(a ! + b t + ..) 
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FIG. III.4: Dynamical diagrams at 2-loop order of type A with two 
external response fields (two connected components) correcting the 
disorder; derived from the two replica static diagrams of Fig. III. 3. 

= X^ a * +bj + •■) ' 

where the combinatorial factors are: 1/3! from the Taylor- 
expansion of the exponential function, 2/2 3 from the explicit 
factors of 1/2 in the interaction, a factor of 3 to chose the 
vertex at the top of the hat, and a factor of 2 for the possible 
two choices in each of the vertices. Furthermore below some 
additional combinatorial factors are given: a factor of 2 for 
generic graphs and 1 if it has the mirror symmetry with respect 
to the vertical axis: each diagram symbol (a.j..) denotes the 
diagram including the symmetry factor. 

We recall that we have defined saturated vertices as vertices 
evaluated at u = while unsaturated vertices still contain u 
explicitly. Diagrams with response-functions added to unsat- 
urated vertices can be obtained by deriving static diagrams: 

a 2 + a 3 = second derivative of the statics 
bi + b 2 + b 3 + b 4 + b 5 + b 6 

= second derivative of the statics . (III.5) 

The graphs which contain external response-fields on satu- 
rated vertices cannot be derivatives from static ones. For class 
A, the hat-diagrams, the only non-zero such graph is ci. 

Explicitly, this reads 

a 2 + a 3 = -d 2 u [-i?"(0)i?"» 2 ] I A , (III.6) 
where (see (A. 18)) 

f d%_d%_ 1_ _L 1 

A : J {2v) d (2n) d q\ + m?ql+ to 2 ((qi+q 2 ) 2 + to 2 ) 2 
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2e 2 



4e 



0(e 2 ) {ehf . 



(III.7) 



Furthermore, we find 



b 4 = [i?»i?'» 2 ] I A (III.8) 



i=l 



and 



ci = 2A'(0+r A"(u)/ A . 



(HI.9) 



The diagram c x is an explicit example for the appearance of 
non-trivial sign-functions resulting from the monotonic in- 
crease of the displacement. It was already discussed in the 
previous Section. In the quasi-static depinning limit (11.11) 
gives (details are given in appendix A): 

Cl = A'(0+) 2 A"(u) 

*/ / 

J qi,Q2 - / *l,*2,*3,*4>0 

e -l(q 2 1 +m 2 )t 1 +{q*+rn 2 )t2 + ((qi+q2) 2 +m 2 )(t 3 +t 4 ) 

sgn(i! - t 3 + U) sgn(t 2 -U + t 3 ) . (III. 10) 
The result of the explicit integration is: 



ei = A'(0+) 2 A'» [J, - I A + finite] 



(III. 11) 



qug2 (<?i + m 2 )(q 2 2 + m 2 )(ql + m 2 )(g 2 + q\ + 2m 2 ) 

(III. 12) 



In 2, x9 
— (e/ 1 ) 2 + finite. 



The last diagram f 2 also involves a sign-function and reads: 

f 2 = 2A'(0+) 2 A'» / / sgn ( t4 -t 3 -t 2 ) x 

•'9li l ?2 Jt 1 ,t 2 ,t 3 ,t 4 >0 
xe -l{qi+q2) 2 +m 2 )(t 3 +t i ) + (q 2 +m 2 )t 1 + (q 2 +rn 2 )t 2 ] 

= -A'(0+) 2 A'»/ ; . (IH.13) 

In appendix A we show that (for any given elasticity) the sum 
of Ci + f2 only involves the integral I a, and that the combina- 
tion takes the simpler form 



ei + f 2 = -A'(0 + ) 2 A»/ A 



(III. 14) 



We now turn to graphs of type B (bubble-diagrams). 
Again diagrams, which are odd functions of u vanish. 
These are: 

hi = h 2 = ii = ji = k 2 = k 3 = 1 2 = 1 3 = U = . (111.15) 
Two other diagram mutually cancel: 

ki+li=0, (IH.16) 
as discussed in Appendix K. 




FIG. III. 5: 2-loop dynamical diagrams of type B (see figure III. 2). 

The diagrams that are second derivative of the static have 
all their response-fields on their unsaturated vertices. These 
are: 



gl + g2 + g3 + g4 = <9 2 



A(u) 2 A'» 



h 3 + h 4 + h 5 + h 6 = d 2 u -A(0)A(u)A"(u) 



12 = ,] 2 



-A(0) 2 A"( U ) 



The surprise is that i 3 , which is not the second derivative of 
a static diagram, since it has both response fields on saturated 
vertices, is non-trivial: 



i 3 = -A'(0+) 2 A»/ 2 . 



(III. 17) 



To summarize, for the driven problem at T = in perturba- 
tion of A = A(m), the contributions to the disorder to one and 
two loops, i.e. the corresponding terms in the effective action 
T[u, u] are: 



^A(u) 
<5 2 A(w) 



A'{u) + (A(«) - A(0))A A 

2l 



(III. 18) 



(A(«)-A(0))A'(«) 
1 



Ia 



J(A( U )-A(0)) 2 A"( U )] I\ 
+A'(0+) 2 A"( u )(/ A -/ 2 ). 



(in. 19) 



Curiously, even though two diagrams contain contributions 
proportional to /; ~ In 2, these contributions cancel in the 
final result for the corrections to the disorder. 

B. Corrections to the friction -q 

We now calculate the divergent corrections to rj, which will 
require a counter-term proportional to iuii. Let us illustrate 
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FIG. III.6: 1-loop dynamical diagram correcting the friction. 



their calculation at leading order. We start from the first order 
expansion of the interaction e -Sin \ which can be written as 



L 



iu xt A(u xt - u xt >) iu xt > ■ 



(111.20) 



>t>t',x 

Contracting one iu xt > leads to: 



/ iu xt A'(u xt - u xt >)R r=0 , t _t> . (IH.21) 

Jt>t',x 

The response function contains a short-time divergence, 
which we deal with in an operator product expansion. Ex- 
panding A' (u xt — u x t> ) to the necessary order yields 

f iu xt [A'(0+) + (u xt -u xt ,)A"(0+) + ...] i? r=0 ,t_ t < . 

Jt>t'.x 

(in.22) 

The first term of this expansion, proportional to A'(0 + ), is 
strongly UV-divergent and non-universal and gives the criti- 
cal force to lowest order in disorder. Since we tune / to be 
exactly at the depinning threshold we do not need to consider 
it. The second contribution, proportional to A"(0 + ), corrects 
the friction: due to the short-range singularity in the response- 
function, we can expand (u xt —u xt ') in a Taylor-series, of 
which only the first term contributes. (III. 22) becomes: 

( iu xt [{t - t')u xt + 0(t - t'f] A"{0+)R r =o,t-t> ■ 

Jt>t',x 

(IH.23) 

The correction to friction at leading order thus is 

Sfi = -A"(0+) Jm r=0}t . (HI.24) 

Here, the response-function is taken at spatial argument 0. In 
momentum representation, the same expression reads 



6 V = -A"(0+) J J tR q . t 
= -A"(0+) f tc^ 2+ 

J a 



m 2 ) 



-A"(0+) 



-A"(0 + )/i 



(q 2 + m 2 ) 2 



(111.25) 



with the already known integral Ji, equation (III. 1). 

We now turn to the 2-loop corrections. There are seven 
contributions, drawn on figure III. 7. Their contribution to rj is 

577 = -ix4x2[a + b + c + d + c + f+g] . (111.26) 



2 for the symmetry of diagrams a, b, e, f and g. Details of 
the calculation of diagrams a to g are given in appendix D. 
Grouping diagrams, which partially cancel, we find: 



a + g=-A"(0+) 2 / 2 



b + c + d = -iA'"(0+)A'(0+)/ 1 2 



(IH.27) 
(111.28) 

o = -A"'(0 + )A'(0 + )J„ (111.29) 
f = -2A'"(0 + )A'(0 + )/ A - 2A"(0+) 2 /4 (IH.30.) 
This involves the non-trivial diagram I v 

1 

(q 2 + m 2 ){q 2 + m 2 ) 2 (q 2 + q 2 + 2m 2 ) 



1 1 - 2 In 2 



2e 2 4e 
calculated in appendix E. 



(eh) 2 + finite (IH.31) 



C. Renormalization program to two loops and 
calculation of counter-terms 

1 . Renormalization of disorder 

Let us now discuss the strategy to renormalize the present the- 
ory where the interaction is not a single coupling-constant, but 
a whole function, the disorder-correlator A(u). We denote by 
Ao the bare disorder - this is the object in which perturba- 
tion theory is carried out - i.e. one consider the bare action 
(II.2) with A — > Ao. We denote here by A the renormal- 
ized dimensionless disorder i.e. the corresponding term in the 
effective action T[u, u] is m e A. 

We define the dimensionless bilinear 1-loop and trilinear 2- 
loop symmetric functions (see (III. 18) and (III. 19)) such that: 



5 (1) (A,A) = m e S 1 A 
S^(A,A,A) = m e S 2 A 



(11132) 
(111.33) 



thus extended to non-equal argument using f(x,y) := 
5 [f{x + y,x + y)- f(x, x) - f(y, y)] and a similar expres- 
sion for the trilinear function. Whenever possible we will use 






, 1 




The combinatorial factor is 1/8 from the interaction, 4 from 
the time-ordering of the vertices, and an additional factor of 



FIG. III.7: 2-loop dynamical diagrams correcting the friction. They 
all have multiplicity 8 except (c) and (d) which have multiplicity 4. 
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the shorthand notation tfW (A) = <5 (1) (A, A) and 6<- 2 \A) = 
S^ 2 '(A, A, A). The expression of A obtained perturbatively 
in powers of A at 2-loop order reads: 

A = m" e A + (5 (1) (m- e A ) + (5 (2) (™- e A ) + 0(A*) . 

(in.34) 

It contains terms of order 1/e and 1/e 2 . This is sufficient to 
calculate the RG-functions at this order (In principle, one has 
to keep the finite part of the 1-loop terms, but we will work 
in a scheme, where these terms are exactly 0, by normalizing 
all diagrams by the 1-loop diagram). Inverting this formula 
yields: 

A = m £ [A - 6 (1) (A) - <5 (2) (A) + ^(A) + . . .1 , 

(in.35) 

where c^ 1,1 ) (A) is the 1-loop repeated counter-term: 



S (hl) {A)^2S^(A,S^(A,A)) 



(111.36) 



The /3-function is by definition the derivative of A at fixed A . 
It reads: 



-md m A\ Ao = e 



m- t Ao + 26 (1) (m- e A ) 



+3(5^ (m- £ A ) + 



(111.37) 



Using the inversion formula (111.35), the /3-function can be 
written in terms of the renormalized disorder A: 



-md m A\ An = e 



A + £«(A) 

+2SW{A)-S ( - 1 ' 1 '>(A) + ... 



(111.38) 



In order to proceed, let us calculate the repeated 1-loop 
counter-term <5 1,1 (A). We start from the 1-loop counter-term 
(III. 18), which has the bilinear form 



2f(u)g'(u) + (/(«) - f(0))g"(u) 

+ (g(u)-g(0))f"(u)]h 



with the dimensionless integral 1\ := l\ 



(IH.39) 
we will use the 



same convention for I a ■= Thus J 1,1 (A) reads 

<5 (1 - 1) (A(u)) = 25« (a,<5«(A)) 

= '(A(u)-A(0)) 2 A'» 

+(A'H 2 -A'(0) 2 )(AH-A(0))]"/ 2 

(in.40) 

Note that this counter-term is non-ambiguos for u — > 0. Fi- 
nally, as discussed at the end of the previous section at any 
point we can rescale the fields ubym^. This amounts to write 
the /3-function for the function A(u) = to~ 2< » A(um'») which 
will be implicit in the following (in addition we will drop the 
tilde superscript). 



The 2-loop /3-function (III. 38) then becomes with the help 
of (111.40) 

-md m A(u) = (e - 2£)A(u) + (uA'{u) 

-i[(A( u )-A(0)) 2 ]" (eh) 

+ [(A( u ) - A(0))A» 2 ]" e(2I A - if) 

+A'{Q+) 2 A"{u)e(2I A -ll^ . (111.41) 

One of our main results is now apparent: the l/e 2 -terms can- 
cel in the corrections to disorder. If it had not been the case it 
would lead to a term of order 1/e in the /3-function and thus 
to non-renormalizability. Thus the /3-function is finite to two 
loops a hallmark of a renormalizable theory. Note that this 
happened in a rather non-trivial way since it required a consis- 
tent evaluation of all anomalous non-analytic diagrams. Fur- 
thermore the precise type of cancellation is unusual: usually 
the 2-loop bubble diagrams of type B are simply the square of 
the 1-loop ones. Here the easily missed and non-trivial bubble 
diagram i 3 was crucial in achieving the above cancellation. 

In order to simplify notations and further calculations, we 
absorb a factor of el\ in the definition of the renormalized 
disorder (or equivalently in the normalization of momentum 
or space integrals). With this, the /3-function takes the simple 
form: 

-md m A(u) = (e - 2QA(u) + C«A'(u) 
-1[(AH-A(0)) 2 ]" 

+ i[(AH-A(0))A» 2 ]" 

+ Ia'(0+) 2 A'» . (IH.42) 

Note several interesting features of this 2-loop /3-function. 
First it contains a non-trivial so called "anomalous term" (the 
last one) which is absent in an analytic theory. Second, it 
can be shown to exhibit irreversibility, precisely due to this 
term. Although, surprisingly, it can be formally be integrated 
twice over u the resulting flow equation for the double prim- 
itive of A(u) does not, however, have the required property 
for the flow of a potential function, i.e. a second cumulant 
of the random potential in the static. This will be shown in 
details in the next Section IV where we find that the fixed 
points of the above equation are manifestly non-potential. In 
Ref. 33 we have obtained the corresponding beta function for 
R(u) in the statics. The corresponding force force correlator 
A s tat(w) = —R"(u) obeys the same equation than (111.42) 
but with the opposite sign for the anomalous term! This shows 
that statics and depinning are indeed two different theories at 
two loops. 

2. Renormalization of friction and dynamical exponent z 

In section IIIB, we have calculated the effective (renormal- 
ized) friction coefficient r)R as a function of the bare one 770 
and the bare disorder A : 



(111.43) 
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This identifies the renormalization group Z-factor as 

Z-^m^Ao] - 1 - Ao(0 + )/i 

+ [A^(0+)] 2 [/ 1 2 + 2/ A ] 
~1 



+ A™(0+)A^(0+) 



-If + 21 A + I„ 



The dynamical exponent z is then given by 



z = 2 + to— In Z(m~ e A ) 
dm 



(IH.44) 



(111.45) 



Equation (111.44) yields 

InZ- 1 = - A / '(0+)/i 
+ A^'(0+) 2 



-If + 21 A 



+ A^"(0 + )A' (0+) 



-I\ + 21 A + Ir, 



(111.46) 



and thus (remind that I x ~ m e and Ia ^ I n ^ m 2e ) 

m^lnZ- 1 = A^O+Xe/x) - A^'(0+) 2 e(/ 1 2 + AI A ) 
+ A^"(0+)A^(0+) e(7 2 + 41 a + 2J„) 

(ni.47) 

We now have to express A in terms of the renormalized dis- 
order A using (III. 35). For the second-order terms, this rela- 
tion is simply Ao = to c A. The non-trivial term is A"(0 + ). 
Using (III. 18), derived twice at + , we get (with the factor of 
(eli ) absorbed into the renormalized disorder 

A^(0+) = (6/ 1 )- 1 x 

(A"(0 + ) + i r i(4A'"(0+)A'(0+) + 3A"(0+) 2 ) 

(IH.48) 



Putting everything together, the result is 



dm 



+e 



4/a 



A"(0+) 2 



2I„ 



?-(S?y-(sf)J A '" <0+ » A ' (0+ > 

(in.49) 



Again there is a non-trivial cancellation of the 1/e terms, an- 
other manifestation of the renormalizability of the theory. In- 
serting the values of the integrals 1a and the dynamical 
exponent z becomes: 



z = 2- A"(0+) + A"(0+) 2 + A'"(0+) A'(0+) 



3 



In 2 



(111.50) 



D. Finiteness and scaling form of correlations and 
response functions 

To complete the 2-loop renormalizability program one must 
check that all correlation and response functions are rendered 
finite by the above counter- terms. In a more conventional the- 
ory that would be more or less automatic. Here however there 
are additional subtleties. The disorder counter-term is a full 
function and is purely static. This counter-term, and its asso- 
ciated FRG equation (111.42) cannot be read at u = because 
of the non-analytic action (this point is further explained in 
Appendix K). Indeed, this equation and the cancellation of 
divergent parts was established only for u ^ 0. It remains to 
be checked that irreducible vertex functions which are u = 
quantities are also rendered finite by the above static « / 
counter-terms. 

We first examine the two point correlation function. We will 
first show that it is purely static. Then, in appendix K we show 
that it is finite and perform its calculation in the renormalized 
theory. One has 



(u 



quj 1 



— V V 



(111.51) 



where lZ quJ is the (exact) response function. We will thus only 
compute T i ui{ l (qt) (in time variable). The 1-loop counter- 
term for rj is absent in this C(A 2 ) calculation of the proper 
vertex but it enters the calculation of (u qix] u^ q - u ) (it dresses 
the external legs R qix] into TZ quJ ). In fact since we find that 
Tiuiu(qt) is static (independent of t) we will need only the 
exact response at zero frequency, which is the bare one be- 
cause of STS. 

To one loop, the proper vertex Tiuiu(qu>) is the sum of the 
graphs a, b, c and d of Fig. III. 6 evaluated at finite frequency 
and momentum, so we write Tiuiu(quj) = a + b + c + d. 
The sum a + b yields after two Wick contractions and short 
distance expansion a term proportional to 



/ 

Jk,ti 



iu t iu t ,A"(u t - u t ')A(u tl - u t2 ) 

*{R k ,t>-t 2 - R k ,t-t 2 ){Rk,t- tl - Rk,t'- tl ) , (HI.52) 

where we have kept all times explicitly to resolve any ambigu- 
ity. Expressing A in a series as in (II. 9), the lowest order term 
is purely static (since one can integrate freely over t\ , t 2 ), and 
proportional to A"(0+)A(0) f k k~ 4 , but vanishes from the 
cancellation between graphs a and b. As explained in detail in 
Appendix K there can be a priori another contribution coming 
from 2A'(0 + ) 2 (5(u)w in the expansion of A"A. It produces 
a term 5(v(t — t'))v\ti — t 2 \ which vanishes when multiplied 
with the above response function combination (since it van- 
ishes at t — t'). 

Thus the only contribution comes from c + d. There the A' 
yields sign functions and there are no ambiguities. One finds: 



d = 



-2A'(0+) 2 f 

J T 



[sgn(i - t 2 ) + sgn(-t - r 2 )] 



1"1,T2>0 



x / e k T ^ c - {k+q) Tl 



=- 2A ' (o+) 7 fe ^ 



q) 2 



-k 2 \t\/ V 
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b = A'(0+) 2 f sgn(r 1 + t/r,) sgn(r 2 - t/rf) 

x / e- fc2T2 e-( fc +«) 2ri 
Jk 

where we have accounted for the extra combinatoric factor of 
2 for graph d and used 

/ c-^sgn(r - t) = \ (6(t)(2c-i 2t - 1) + 0(-tj) . 
Jt>o Q v 7 

(in.53) 

We thus find that although each graph is time dependent, 
this time-dependence cancels in the sum. Thus we find a static 
result: 



A(0) 



A'(0+) 2 f ■ 



l 



(IH.54) 

The static 1-loop counter-term should thus be sufficient to 
cancel the divergence of (111.54). This is further analyzed in 
Appendix J where the full correlation-function is computed. 

We have thus found the commutation Tim^u = Q,q) = 
Tiuiu(u = + ,q). Note that if all correlation functions are 
purely static, i.e. strictly time-independent, it implies the com- 
mutation of the limits. Then it also implies the finiteness 
since these static divergences have been removed. We have 
not pushed the analysis further but we found a simple argu- 
ment which indicates that all correlations are indeed static. 
We found that the time dependence in diagrams cancels by 
subsets, noting 63 that graphs can be grouped in subsets (e.g. 
pairs ac, bd, cf in Fig. III. 3) which vanish by shifting the end- 
point of an internal line within a splitted vertex. 

Finally, let us note that our result that correlations at the 
quasi static depinning are purely static for v = + is at vari- 
ance with previous works 19 ' 20 . Thus the only functions where 
the dynamical exponent comes in are response function. 

E. Long range elasticity 

As was discussed in the Introduction there are physical sys- 
tems where the elastic energy does not scale with the square 
of the wave-vector q as Elastic ~ q 2 but as -Elastic ~ q a - In 
this situation, the upper critical dimension is d c = 2a and we 
define: 



e := 2a - d . 



(111.55) 



The most interesting case, a priori relevant to model wetting 
or crack-front propagation is a = 1, thus d c = 2. 

In order to proceed, we have again to specify a cut-off pro- 
cedure. For calculational convenience, we choose the elastic 
energy to be 



E, 



clastic 



(q 2 + m 2 )% . 



This changes the response-function to 

R qt = e(t)e^ 2+m2 ) ft 



(111.56) 



(IH.57) 



Since contributions proportional to see (A. 26), cancel, the 
only integrals which appear in the /3-function are: 



1 J (9 2 +'. 



9 ( 9 ^r m wi eV (IIL58) 



1 



7 («) = 

'51,92 {qi + m 2 )^{ql+m 2 ) a {{qi+q2) 2 + m 2 ) L : 

(IH.59) 

The important combination is again 21 { ^-{l[ a) ) 2 . One finds 
(see appendix F) 



2e{2I^~{l[ a) f) 
1 dt 1 + -(l + t)j_ + T'{a) r'(f) 



Jo ~t 



lo t (1 + t) 
0(e). 



T(a) r(f) 

(in.60) 



Since this term is finite, the /3-function is finite; this is of 
course necessary for the theory to be renormalizable. For the 
cases of interest a = 1 and a = 2, we find 

X {2) = 1 (IH.61) 
X (1) =41n2. (IH.62) 

Since there is only one non-trivial diagram at second order, all 
2-loop terms in the /3-function get multiplied by 

-md m A(u) = (e - 2()A(u) + (uA'(u) 
-i[(AH-A(0)) 2 ]" 

+^[(A( U )-A(0))A'( U ) 2 ]" 



+ 



^A'(0+) 2 A"(u) . (111.63) 



The diagrams involved in the dynamics also change. Besides 

1 

5 (g|+m 2 ) [(q|+m 2 )3 + (ql+m 2 )^] 

(ni.64) 



and 1^ given above we need 



/(i) •= f 

" ' J (q 2 +m 2 )h 



1 In 2 - -j \ , m . 2 
— + (el[y+ finite 

calculated in appendix G. 

Starting from (111.49), the dynamical exponent z is then in 
straightforward generalization of (III. 50) given by 

z = a - A"(0+) + A^A"(0+) 2 + y( a >A"'(0 + )A'(0 + ) 

(IH.65) 

with given above and 



y(«) — + 



24 a) - (l[ a) 



r( 1 )=6In2-- 
2 

y( 2 ) = --ln2. 
2 



(in.66) 

(111.67) 
(111.68) 
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The case a — 2 reproduces (III. 50). Since both and 

are finite, we have checked that also in the case of long-range 

elasticity the theory is renormalizable at second order. 

IV. ANALYSIS OF FIXED POINTS AND 
PHYSICAL RESULTS 

The FRG-equation derived above describes several different 
physical situations: periodic systems (such as charge den- 
sity waves) where the disorder correlator is periodic and non- 
periodic systems (such as a a domain-wall in a magnet). 
Within the latter, SR (random bond) and LR (random field) 
disorder must a priori be distinguished. In our analysis of the 
FRG-equations, we have to study these situations separately. 

Before we do so, let us mention an important property, valid 
under all conditions: If A(u) is solution of (III. 63), then 



A(u) := k 2 A(u/k) 



(IV. 1) 



is also a solution. We can use this property to fix A(0) in 
the case of non-periodic disorder. (For periodic disorder the 
solution is unique, since the period is fixed.) 

A. Non-periodic systems 

We now start our analysis with non-periodic systems, either 
with random field disorder or any correlator decreasing faster 
than RF. Let us first recall that at the level of the bare model 
the static RF obeys R(u) ~ — a\u\ at large \u\ and thus 
J °° duA(u) = R'(0+) - R'(oo) = -a (a is the amplitude of 
the random field) while RB or any correlator decaying faster 
than RF satisfies / A = 0. 

Let us first integrate the disorder flow equation (III. 63) from 
u = + to u = +oo. We obtain 



-md. 



oo oo 

. J A{u) du=(e- 3C) J A(u) du - X^A'(0+ f . 



(IV.2) 

The only assumption that we have made here is that uA(u) 
goes to zero at u = +oo, which is the case both for RB and 
RF. 

Let us first recall the 1-loop analysis, where in the FRG 
equation there is no distinction between statics and depinning. 
The last term in (IV.2) is then absent. Thus one finds either 
fixed points with J A = a > with £ = e/3, the RF univer- 
sality class, or others with J A — for £ < e/3 which cor- 
responds to disorder with shorter range correlations than than 
RF. This includes the RB fixed point with exponentially de- 
caying correlator and Qrb = 0.208e. It also includes a contin- 
uous family of intermediate power law fixed points 32 ' 36 with 
decay at large u as A(u) = —R"(u) <~ (a — 2)(a — l)u~ a 
with a* > a > 1. These have ((a) — e/(2+a) (from solving 
the linear part of the FRG equation) and ((a*) = (re- 

The last term in (IV.2) shows that things work differently to 
two loops at depinning. the condition J A = is no longer 
possible at the fixed point. Starting from RB one develops a 
positive value for J A, i.e. a random field component. The 



natural conclusion is then that all correlations shorter range 
than RF flow to the RF universality class 76 . Furthermore the 
fixed point condition (III. 63) equals implies a unique well 
defined value for £ identical for all ranges shorter than RF 
(including RB). This value takes the form: 

C=|- t ( } = \ + ^ + 0(e 3 ) , (IV.3) 

3/A 

o 

where C2 > can be obtained from the knowledge of the 1- 
loop fixed point A ~ 0(e) only. 

Before we compute £ 2 and obtain the depinning fixed point 
to two loop, let us note that in the static case 33 the last term in 
(III. 63) has the opposite sign and, integrating over u one finds 
that there is thus no term proportional to A'(0 + ) in (IV.2). 
Thus for the RF disorder < / A < +00 one can again 
conclude that 



Alt* _ _ 



(IV.4) 



to (at least) second order. In fact, as discussed in 33 this is 
expected to be exact to all orders due to the potentiality re- 
quirement of the static FRG equation, which also implies that 
J A = holds to all orders at the static RB fixed point. The 
corresponding value for is given to order e 2 in 33 . 

We now want to find the fixed-point function of equation 
(111.63). Using the reparametrization-invariance (IV.l), we set 
(with the factors 1/3 and 1/18 chosen for later convenience) 



A(u) 

2/1 (0) 
2/2(0) 



2/1 0) 
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y 2 (u)+0(e 3 ) 



(IV.5) 

(IV.6) 
(IV.7) 



where y\ (u) is the 1-loop fixed point function for the RF case. 
It was obtained in Ref. 22 and further studied in Ref. 28 . Let us 
recall its properties. To lowest order in e the 1 -loop /3-function 
(IH.63) reads: 



-A(u) + -uA'{u) - - 



(A(u)-A(O)) 
inserting (IV.5) the function yi(u) must satisfy: 



, 1 r 

[uyi(u)} = - (yi(u) - 2/1(0))' 



which can be first integrated to 

uyi(u) = (2/1 (u) - 1)2/1 (w) , 



: (IV.8) 



(IV.9) 



(IV. 10) 



using (IV.6) in the last line. A second integration with the 
boundary conditions implied by (IV.6) yields: 



2/i (u) — In 2/1 (u) = 1 + 



1 



(IV. 11) 



The derivatives of 2/1 (u) at u = will be needed below. De- 
riving (IV.10) successively w.r.t. u, we find 



2/i(0) = 1 , 
2/i'(0 + ) - \ , 



2/i(0 + ) = 
2/i"(0 + ) 



-1 



(IV. 12) 
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FIG. IV. 1: The fixed-point function yi(u) at 1-loop order for non- 
periodic disorder. 



We have also determined the fixed-point function at second 
order yi (u), which is given in appendix H. 
In order to extract ( from (IV.2), we need 



/•OO 

V27 = / yi ( 
Jo 



u) du . 



(IV. 13) 



which was computed in Ref. 28 . The method is to convert 
(IV. 13) into an integral over y := y\\ 



dy 



v^Vy-iny-i ' 



(IV. 14) 



where in the last equality (IV. 11) has been used. Integrating 
by parts, this yields 



7 



Jo 



lny-lw 0.5482228893... (IV. 15) 



Combining the definition of A(u) in (IV.5) with (IV. 12) and 
(IV. 15), we find 



(2 



x(<*) 

27V2 7 



and thus for ( 



3 + 27V2 7 



+ 0(6 3 ). 



(IV. 16) 



(IV. 17) 



This result violates the conjecture of , that £ = | to all orders 
in e. To compare (IV. 17) with simulations, we have to specify 
to the cases of interest: First, for short-range elasticity, i.e. 
a = 2, we find 

C = I (l + ^-)=|(l + 0.143313 e). (IV.18) 

Our results are in excellent agreement with the numerical sim- 
ulations, see figure IV.2. For long-range elasticity, i.e. a = 1, 



exponent 


dim 


1-loop 


2-loop 


estimate 
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FIG. IV.2: Depinning exponents, a = 2. First column: Exponents 
obtained by setting e = 1 in the 1-loop result. Second column: Ex- 
ponents obtained by setting e = 1 in the two-loop result. Third col- 
umn: Conservative estimates based on three Pade estimates, scaling 
relations and common sense. Fourth column: Results of numerical 
simulations obtained directly without using scaling relations. 



equation (IV. 17) reads 



c 



4 In 2 

9V2 7 



e^j = 0.39735 e) . (IV. 19) 



This is in reasonable agreement with simulations, as shown on 
figure IV. 3. 

We now turn to the calculation of the dynamical expo- 
nent z. As can be seen from the general result of equation 
(IH.65), we need A'(0+), A"(0+) and A'"(0+) at leading or- 
der, which can be inferred from (IV.5) and (IV.12). We further 
need A"(0) at second order. Expanding (III. 63) to order e 3 , 
and Taylor-expanding to second order in u, we can solve for 
y'2 (0), which yields 



A"(0) = |i,i'(0) + ^(0) 



(IV.20) 



exponents 


one loop 


two loop 


estimate 


simulation 
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0.33 


0.47 
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FIG. IV.3: Depinning exponents, a = 1. First column: Exponents 
obtained by setting e = 1 in the 1-loop result. Second column: Expo- 
nents obtained by setting e = 1 in the two-loop result. Third column: 
Conservative estimates based on three Pade estimates together with 
scaling relations between exponents. 
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(IV.21) 



Specializing to the SR-case (a = 2) yields with the help of 

(in.50) 

z = 2- A"(0+) + A"(0+) 2 + A"'(0+)A'(0+) Q - In 2^ 



2-* 
9 



In 2 



5 

108 



2 - 0.222222e - 0.0432087e 2 



(IV.22) 



The agreement with the numerical simulations given on figure 
IV.2 is again good. Finally, the exponents (3 and v are obtained 
from scaling relations. For a = 2 (SR) they read 



c 



2-C 



2 f C2 

6 



1 - 



9 

1 _ 1 

~C ~ 2 
e 



0.040123e 



1 

24 
2 



ln2\ 
108j 



e 

12 



1 

72 



- + — + 0.0258316 e 2 



(IV.23) 



(IV.24) 



We now turn to long range elasticity a 7^ 2. The general 
formula for z reads 



z = a - -e + e 



3 



2X^ 

+ 



27 



54 



(IV.25) 



Specifying to a = 1 yields 



>4 



7r + 201n2 
108 



= 1 - -e - 0.1132997e 2 . (IV.26) 
Again (3 and v are obtained from scaling as (a = 1) 

z-C, , 1 . o.( Co. 2 tt + 20 In 2 



/3 



i-C 

2 



9 



2 
27 



108 



e- 0.1873737e 2 



Z/= 1 + -+6 2 ( -+C2 

3 V 9 s 



1 



0.24356e 2 



(IV.27) 



(IV.28) 



Numerical values are given on figure IV. 3. 

Note that to two loops at the RF fixed point there does not 
appear to be any unstable direction. We thus conclude, as in 21 
that 



v . 



(IV.29) 



Finally, for depinning there should also be a family of fixed 
points corresponding to correlations of the force which are 



long range with A(u) ~ nT a and a > a* > 1. The linear 
part of the FRG equation implies that ((a) = e/(2 + a) and 
a crossover from the RF fixed point occurs when ((a*) — 
Crf = § + C2£ 2 + 0(e 3 ). We have not studied these LR fixed 
points in details. 

B. Periodic systems 
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2-loop 
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0.83 
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0.78±0.03 


0.81±0.03 
0.84±0.05 
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0.52±0.08 


0.63±0.06 
0.68±0.07 
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FIG. IV.4: Depinning exponents for CDW. First column: Exponents 
obtained by setting e = 1 in the 1-loop result. Second column: Ex- 
ponents obtained by setting e = 1 in the two-loop result. Third col- 
umn: Conservative estimates based on three Pade estimates, scaling 
relations and common sense. Fourth column: Simulations from 9 . 

For periodic A(u) as e.g. CDW depinning 17,21 , there is an- 
other fixed point of (III. 63). It is sufficient to study the case 
where the period is set to unity, all other cases are easily 
obtained using the reparametrization invariance of equation 
(IV. 1). This means, that no rescaling is possible in that direc- 
tion, and thus the rescaling factor is 



C = 



(IV.30) 



The fixed-point function is then periodic, and can in the inter- 
val [0, 1] be expanded in a Taylor-series in u(l — u). Even- 
more, the ansatz 

A(u) = (one + a 2 e 2 + ...) + (V + b 2 e 2 + . . .) u(l - u) 

(IV.31) 

allows to satisfy the fixed-point equation (III. 63) to order e 2 , 
with coefficients 



A*(u) 



e e 
36 + _ 



108 



6 



9 



u(l-u) . (IV.32) 



In the physically interesting situation of charge density waves, 
the elasticity is short range, i.e. a = 2 and = 1 which 
yields: 

A*(u) = 4; + ^:-(^ + ^)u(l-u). (IV.33) 
u 36 108 \6 9j 1 ' 

This fixed point is manifestly non-potential, i.e. it describes 
a force-force correlation-function, where the forces can not 
be derived from a potential. In a potential environment, the 
integral of the force over one period must vanish, and so must 
the force-force correlation-function. In contrast we find here 



./o 



&uA*{u) 



108 



. (IV.34) 

108 v ' 



Thus to two loop the fixed point correctly accounts for the 
irreversibility in the driven system, which becomes manifest 
beyond the Larkin-length. This was not apparent to one loop. 
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An important feature of the periodic case is that the fixed 
point is unstable. The direction of instability is simply adding 
a constant to A(u) and its eigenvalue is trivial equal to e to two 
loops and presumably to all orders. The full stability analysis 
is performed in appendix I but it can be seen already from: 

1 1 
-md m J A(u) du = e f A ( u ) dw - ( Q >A'(0 + ) 3 
o o 

(IV.35) 

obtained by integration of the 2-loop FRG equation on the 
interval [0 + , l - ]. One sees that J A flows away if it does not 
coincide with its fixed point value (IV. 34). 

Thus the asymptotic flow as the dimensional parameter 
in — ► takes the simple form 

A m (u) = A* (u) + cm~ e (IV.36) 
i 

c = m £ J (A mo - A*) (IV.37) 
o 

i.e. it takes the fixed point form shifted by a growing constant. 
In the statics J du A m = f Q du A* =0 from potentiality ( 
the last term in (IV.35) is absent ) and thus c = 0. At depin- 
ning c is non-zero at 2-loop order (c w — J Q A* > 0) using 
that the bare disorder has zero integral) and this has several 
consequences. First one obtains the static deformations as the 
sum 



(u x - u Q ) 2 = B n i(a:) + B RF (x) 
of a universal logarithmic growth-term 



.Bni(ar) = A d In \x\ 



Ad = vh e+ 



2X( Q ) 



108 



(IV.38) 

(IV. 39) 
(IV.40) 



(the calculation of Ad is presented in Appendix J as an exam- 
ple of an explicit calculation of a correlation function in the 
renormalized theory); and of the contribution of the generated 
"random force" of the Larkin type 



B RF (x) ~ c\x\ 



(IV.41) 



which completely decouples from the other one. This is very 
similar to what was found in other driven systems where a 
random force is generated 50 ' 66 . In particular this implies that 
the true roughness-exponent at depinning is not ( = but 



Cdep 



4-d 



(IV. 42) 



Another consequence is that the the two exponents v and 
v FS are different. We find: 



1 



1 



2-C 
1 



^FS 



Cdep 



(IV.43) 
(IV.44) 



and given the generality of the above argument this should 
holds to all orders. Note then that the CCFFS -bound 48 for 
^fs is saturated. This is very different to the case of interfaces 
(saturation of the bound there would lead to the incorrect re- 
sult C = e/3). 

The dynamical exponent z is 



2 

3 9 



(IV.45) 



Curiously, it does not depend on the diagram l n or equiva- 
lent^ rw. 

In CDW depinning, the best observable quantity is 0. From 
the scaling relation 17 - 19 - 21 = (z - ()/(2 - (), and C = 0, 
we find = z/2 and thus for CDW (a = 2) 



0=1- 



18 ■ 



(IV.46) 



This expansion is however ill-behaved, at least at large e. It 
therefore seems advisable, to use one of the Pade-variants. 
The only one which respects common sense down to d = 1 
and even beyond, is the Pade (0,2), reading 



= 



1 



1 + - 4- — 

1 ^ 6 T 12 



(IV.47) 



Again simulations are in reasonable good agreement with our 
theoretical predictions, as can be seen on table IV.4. Further 
simulations would be welcome. 

V. CONCLUSION 

To conclude we have constructed a consistent field theory of 
isotropic depinning at zero temperature to 2-loop order. While 
the 1-loop flow-equations for statics and driven dynamics are 
identical, our 2-loop equations distinguish these physically 
different situations, yielding different universal predictions for 
both cases. This is an encouraging progress. The non-analytic 
field theory that we have developed here will be discussed in 
companion studies 34 35 for the static theory to two and three 
loops. 

A lot remains to be done and understood. If universality is 
to hold at depinning then a renormalizable theory should exist 
to any number of loops. We have not attempted a proof to all 
orders here, and the mechanism in which the 1 / e-divergences 
cancel is non-trivial. We have however checked the applica- 
bility of formal constructions like the subtraction operator R 
on sample diagrams. This could further be tested in a 3-loop 
calculation. Although short time singularities (S(vt) terms) 
did not appear their role to any order remains to be clarified. 
Next, effects of temperature have not been included here. One 
expects that although at T = the statics and the depinning 
should be two distinct field theories, this distinction becomes 
blurred at finite temperature. How this will work out is not yet 
elucidated. Some efforts in that direction are reported in 77 . 
Similarly it would be quite interesting to understand how to 
describe / = /~, i.e. the approach to the threshold from be- 
low. From the considerations here this appears to be quite 
non-trivial. 
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Extension of the present method to systems with N > 1 
is also far from trivial. The monotonous increase of u t 
does not apply to all components, which leads to complica- 
tions. The large- N limit of the static FRG was solved exactly 
recently 36 and it would be interesting to extend it to the dynam- 
ics. Finally the threshold dynamics of other systems, such as 
random field spin models, which can be described by the FRG, 
is of interest. 

From the point of view of simulations our results together 
with recent more powerful algorithms offer hope that more 
precise comparisons could be made, not only for exponents 
but also for other universal quantities which offer stronger 
tests such as scaling functions, amplitudes or finite size ef- 
fects. The exponent z^fs should be measured independently. 
We encourage further precise numerical studies on both man- 
ifolds and CDW with a comparison to theory in mind. Agree- 
ment between numerics and theory would allow to rule out or 
to accept elastic models for the description of more complex 
experimental situations. 

APPENDIX A: CORRECTIONS TO 
DISORDER: DIAGRAMS OF TYPE A 

In the following, we give explicit expressions for the diagrams 
contributing to the renormalization of disorder. To simplify 
notations, we have introduced q 3 := qi — q 2 and set the mass 
m to zero. The mass-dependence can easily be reconstructed 



by replacing qf by qf 



We start with the diagrams of 



class A, given on figure III.4. For illustration, we show the 
complete calculation of the first non-vanishing diagram &2'. 



T * ! * Jqim Jt 1 ,t 2 ,t 3 ,t 4 >o 
tr i 



92*2-9 3 (*3+*<l) 



xA"'(u T - u a )A'(u T - t2 - Uer-ti) 
xA(« T _ t3 -ti T _ f2 _ (4 ). (A.l) 

For the field u, we have given the time-arguments, but sup- 
pressed the spatial arguments, since the result is taken at con- 
stant background field. We also do not write explicitly the two 
response-fields. The given configuration is for r — a large, 



and we can set u := u T i — u a i for all t' = r ± some tj, 
a' = a ± some ti since the ti are small (due to the exponen- 
tially suppressing factors) compared to the difference of r — a. 
Finally, since A is continuous, A(u T - t3 — u T - t2 -t 4 ) can be 
replaced by A(0). Integrating over all times leads to 



(A.2) 



a 2 = -2A(0)A»A"» f * 

J gi, q 2 HlH^ 

Similarly, we find 

a 3 = -2A(0)A» 2 f * 

J 91,92 ^1^2 



a 2 + a 3 = -dl 



A(0)A» 2 



2 4 

i 



2 2 4 
91,92^1^293 



2 2 4 
91,92 



bi = 2A'(u) 2 A"(u) 
b 2 = 2A(u)A'(u)A'» / -A 4 



b 3 = A(u)A"(u) 2 I 

Jqi,q 2 Qlll 



b 4 
b 5 



2 4 

qi,q2 ^1^2% 

A"(u)A'(u) 2 f ' 

Jqi,q 2 I1Q2I3 



2A"(u)A» i 



1 



91,92 9l92?3 



b 6 = A( M )A» 2 f * 

•/«, ,0, 111-; 



^ 2 4 
9l,92 I1I2I3 



The contribution of the bj's can be summed as 

1 



1=1 

Diagram ci is 



A(u)A'(u)' 



91,92 mhi 



(A3) 

(A.4) 
(A.5) 

(A.6) 

(A.7) 

(A.8) 

(A.9) 

(A. 10) 



(A. 11) 



ci = 2A'(0+) 2 A» f (A.12) 

•^91,92 1ll2l3 



All these diagrams contain the hat-diagram known from the statics and </> 4 theory. It can be calculated as follows: 



A l quq2 {ql 



+ m?){ql + m 2 ) 2 (( qi + q 2 ) 2 + m 2 ) J a ^ J>0 



L 



0e 



i q ^ i a 



>0,/3>0,7>0 



f3e 



-m (a+/3+7) 



Det 



a + 7 7 
7 /3 + 7 



-d/2 



3-d -m 2 7 (a+/3+l 



a>0,/3>0,7>0 



\a + l3 + al3)- d ' 2 



e~ q r(4-d)m~ 2e J, 



(A. 13) 
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where we split the divergent integral J in pieces, which are either finite or where the divergence can be calculated analytically: 

(a + + a0) ' 



/>oo />oo 

/ d «/ 

Jo Jo 

POO rl 

/ da 
Jo Jo 



J 
Jl 
J 2 
■h 



d0 
d0 



(a + + a0) 2 (a + 0+iy 
(a + + a0)i 



(a + + a0f {a + + iy 

(a + + a0)i 



= Jl + J 2 + J 3 

= In 2 + 0(e) 



1 



OO /'OO 

^Ji \(a + + a0)2 (a + + iy (l + a) 2 -?0 1+ i 

OO rOO 

da / d0 
Ji 



ln2 + (9(e) 



(l + a)»-^ + t = ! + 1 + Q(e) 



(A. 14) 
(A.15) 
(A. 16) 
(A. 17) 



r 



This gives the final result for the hat-diagram 



I A = 



e~ q T(4 - d) m" 2e - + 1 + (3(e 



1 1 

2e2 + Te 



(A. 18) 



26{t 3 - U)e- q2 ^ t3 - u) + 29(U - t 3 )e^ 2 (*4-t 3 ) - 1) 

(A.22) 



The integral over the two remaining times t 3 and t 4 in (A. 20) 
gives 



We now turn to the non-trivial diagram ei . At finite velocity 
v, the diagram is 



/ / A"(u)A'(v(ti+t 4 -t 3 )) 

J 91,92 J tl,...,*4>0 

xA'(v(t 2 + t 3 - ti))e- tiq *- t2q *- {t3+u)q * 

(A. 19) 



/ I(t 3 ,t 4 ,q 1 ,q 2 ) = 

Jt 3 ,t4.>0 . 




Q1Q2 \3s(Si 5f(92+93) 4 



(A.23) 



and thus 



In the limit of vanishing velocity v — > 0, we can replace 6l = A'(0 + ) 2 A"(w) x 
A'(v(t 2 +t 3 - t 4 )) by A(0+)sgn(t 2 + h - h) a.s.o. Let 
us stress that this replacement is correct both before and after 
reaching the Larkin length. Its result is 



x 



/ -A 



1 



1 



1 



'91 ,92 9l9293 W +93 9^+ 93 93 

with (j 3 = qi — q 2 . In presence of a mass this reads: 



ei = A'(0+) 2 A"( 



11 



(A.24) 



< / / e- q ^ t3+t ^I(t 3 ,t 4 , qi ,q 2 ) (A.20) 

'91,92 Jt3,t4>0 

I{t 3 ,t 4 , qi ,q 2 )= f 0(ti)0(t2)sgn(ii+t4-t3) 

"'tl,*2 

x sgn(t 2 +i 3 - t 4 )c _(9 i* 1+922 * 2) . 

(A.21) 

Using that c-(^ti+9 2 2 t2) = a e -(«i' 1+ ^ ! ' and inte- 

& g 2 g 2 Cl I 2 

grating I by parts in ti and t 2 yields 



ei = A'(0+) 2 A'» x 

x 



1 



91,92 (9 2 +™ 2 )(9 2 2 + ™ 2 )( g 2 + m 2 ) 
1 1 



+ ■ 



1 



I(t 3 ,t 4 ,qi 1 q 2 ) = 



1 



qi<& 



q\ + ql + 2m 2 q\ + q\ + 2m 2 g 2 + to 2 / 

(A.25) 



We now calculate the new integral. It is relatively simple, 
since it has only a single pole in 1 /e: 



, /uq2 (qi + m?){ql + m 2 )(g 2 + m 2 )(g 2 + q\ + 2m 2 ) 

/ e -a(q 2 1 + m 2 )-l3((q 1 +q 3 ) 2 + m 2 )-y(qi+m 2 )-8(q 2 1 +q 2 3 +2m 2 ) 
91,93 Ja>0,(J>0,7>0i>0 
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-m 2 (a+/3+7+2<5) 



q>0,/3>0,7>0,<5>0 



Det 



a + p + 5 (3 
(3 /? + <5 + 7 



-d/2 



e - m 2 i(a+/3+ 7 +2) (5 3-^ 1 +1 + a + 2/3 + a[3 + a~f + /3 7 )" 



d/2 



a>0,/3>0,7>0,i5>0 
2 



= r(4-d) 

Mi 

— (e/i) 2 + finite 



(1 + 7 + a + 2/3 + a/3 + aj + /?7)" d/2 [m 2 (a + (3 + 7 + 2)] 



d-4 



= 21n2 



a>0,/3>0,7>0 

(1 + 7 + a + 2(3 + a[3 + 07 + /?7)~ 2 + finite 

a>0,,3>0,7>0 

2e + finite 



(A.26) 



This gives the final result for ei 

ei = A'(0+) 2 A'» (2J ; - 7 A ) . 
The last non-vanishing diagram is f 2 : 

e -9l(*3+*4)-(9?tl+g2*2 



(A.27) 



f 2 = 2A'(0+) 2 A'» [ 
Jgi 



92 



x sgn(t 4 - t 3 - t 2 ) . (A.28) 



Integrating first over t± and then over the remaining times 
gives 



-2A 



'(0 + ) 2 A»/ - :7 ;- :;:: ? - ^ 



(A.29) 



The integral has already been calculated in (A.26), yielding 
the result 



f 2 = -2A'(0+) 2 A'»7, 



(A.30) 



Note that the non-trivial integrals in ei and f 2 are in fact iden- 
tical and cancel: 



e 1+ f 2 = -A'(0 + ) 2 A»/.4. 



(A.31) 



APPENDIX B: CORRECTIONS TO 
DISORDER: DIAGRAMS OF TYPE B 

In this appendix, we calculate diagrams of type B (the bubble- 
chains). 

The diagrams which are odd functions of u are: 



hi 



: Ji 



12 = 13 = 14 = 0. (B.l) 



The diagrams that are second derivative of static ones have all 
their response-fields on their unsaturated vertices. These are: 

gl = A'» 2 A/ 2 (B.2) 
g 2 = 2A»A"»A(u)/ 2 (B.3) 



g 3 = A» 2 A»/ 2 



iA(«) 2 A"»I? 



;i + g2 + g3 + g 4 = a 2 



-A(u) 2 A'» 



h 3 

h-5 = h 6 
h 3 + h 4 + h 5 + h 6 

12 = h 
ki = -li 



(B.4) 
(B.5) 

(B.6) 



-A(0)A"'»A(u)/ 2 (B.7) 
-A(0)A'» 2 / 2 (B.8) 
-A(0)A"»A»J 2 (B.9) 

d 2 [-A(0)A( u )A»]/ 2 (B.10) 

iA(0) 2 A""(u)/ 2 (B.ll) 

-A(u)"A"(0+)A(0)/ 2 . (B.12) 



The surprise is that i 3 , which is not the second derivative of 
a static diagram (since it has both u on saturated vertices) is 
non-trivial: 

i 3 = -A / (0+) 2 A"(u)7 1 2 . (B.13) 
This diagram is necessary to ensure renormalizability. 

APPENDIX C: CORRECTIONS TO 
DISORDER: DIAGRAMS OF TYPE C 

In this appendix, we show that diagrams of type C do not con- 
tribute to the renormalization of disorder. This is fortunate, 
since they involve a strongly diverging diagram (the tadpole), 
which would render perturbation theory non-universal. 
The diagrams which are odd functions of u are 

mi = rn 3 = m 4 = m.5 = ni = n 3 

= n 4 = n 5 = p 2 = p 3 = q 2 = q 3 = . (C. 1) 

The following diagrams cancel: 

m 2 + n 2 = (C.2) 

Pi + qi = (C.3) 

P4 + 94 = . (C.4) 
No contribution remains. 
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ql 



3, 




q2 



FIG. C. 1 : 2-loop diagrams of class C 

APPENDIX D: CORRECTIONS TO rj: 
2-LOOP DIAGRAMS 

In this appendix, we give all diagrams contributing to the cor- 
rection of T] at second order. For simplicity of notation, we 
again drop the explicit mass-dependence. We group together 
those diagrams which partially cancel. We demonstrate ex- 
plicitly how to calculate the first diagram a from the very be- 
ginning. 



(D.l) 



RqitiRq 2 t 2 Rqit 3 X 




91,92 t 1 ,t 2 ,t 3 



A'K=o - u_ tl _ t2 _ t 3)(-A"(ti_ tl )) 



We have drawn three response functions. We have chosen 
to start counting time at for vertex a, such that vertex (3 
is at time —t\, vertex 7 at time —ti — t 2 and vertex S at time 
—ti — t 2 — t 3 . This gives the times for the arguments of A. The 
upper A in diagram a has one time derivative, the lower vertex 
two, resulting in A' and — A" respectively (the minus-sign is 
a consequence of the both response-functions entering at dif- 
ferent "ends" of A). We have suppressed the space-arguments 
in the fields u, since all diagrams correcting 77 are calculated at 
a spatially constant background field. Inserting the response- 
functions R qt — Q{t)e~ q *, we arrive at 



■/ / 

91,92 tl,*2,t3>0 



-qf(ti+t 3 )-qit 2 



A'(u - ti_ tl _ t2 _ t3 )A"(u_ tl - u- tl - t2 ) 

(D.2) 

The crucial point is now that this diagram corrects the critical 
force and r\. The correction to the critical force is obtained by 
setting the arguments of the A's to + (unique here due to the 
time-arguments). This contribution is non-universal and we 
shall not calculate it in the following. The universal correc- 
tion to 77 is obtained by Taylor-expanding the argument of e.g. 
A'(u - w_ tl _ t2 _ t3 ) as 

u - w_ tl _ t2 _ t3 w (v + u )(h + t 2 + t 3 ) (D.3) 

and thus 

A' (m - u_ tl -t 2 -t,) « A"(0 + )uo(ti + t 2 + is) , (D.4) 

which naturally leads to the generation of a correction to fric- 
tion. For our diagram, this is (sloppily dropping u and the 
response-field for simplicity of notation) 



a = 



_ J J e -ql(ti+t 3 )-qit 2 x 

91,92 tl,*2,*3>>0 

[A"(0+) 2 (t! + h + t s ) + A'(0+)A"'(0+)t 2 ] 



J A"(0+) 2 ( 



r + - 



fi 2 1 4 4 

qm Q1I2 



+ A'(0+)A"'(0+) 



4 4 



91,92 



(D.5) 



where in the last line we have explicitly performed the time- 
integrations. 
Diagram g is 



-9?(tl+«3)-92*2 



91,92 ti,t 2 ,t 3 >0 



A'(u - u- tl - t3 )A"(u- tl - u- tl - t2 ) 

= J j e -qf(ti+t 3 )-qlt 2 x 

91,92 tl,*2,*3>0 

(A"(0+) 2 (ti + i 3 ) + A'(0+)A"'(0+)i 2 ) 
= / A "(° + ) 2 i + A'(0+)A"'(0+)-l I . (D.6) 

Jqi,q 2 ^1^2 ^2 

Thus 

a + g=- I A"(0 + ) 2 ^4 = -A"(0+) 2 / 2 (D.7) 

J 1\l2 

91,92 

Note that both diagrams a and g contain a tadpole-like sub- 
divergence, which is canceled by a counter-term for the crit- 
ical force. However their sum does not involve such a term 
and thus there is no need to specify it. 
Graphs b, c and d: 

b = _ y j e -9?(ti+*3)-92*2 x 

91,92 tl,*2,*3>0 
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A'"(m - u_ t2 )A(u_t 3 - u- tl - 



12) 



e -qi(ti+t 3 )-q 2 t 2 x 



91,92 J ti,t 2 ,t 3 

mi /r\+\ A /'n^-^ 1 A "' 



A""(0+)A(0)t 2 + A"'(0+)A'(0 + )|ti +t 2 -ti 



(D.8) 



while 



+ «3)-<?2*2 



91,92 *l,t 2 ,*3>0 

A""(0+)A(0)t 2 + A"'(0+)A'(0+)|t 1 - i 3 | 



Note the factor 1/2 for the symmetry in the lower vertex. To- 
gether they are 

b+c + d = A'"(0+)A'(0+) 



91,92 ti,t 2 ,t 3 



(D.10) 



Changing variables to u = t\ — i 3 and s = 5(^1 + £3), the 
integral J ti t ^ >Q becomes J °° ds J_ s 2g du. The integral over 



u can be performed, but for fixed s the second term in (D. 10) 
depends on the value of t 2 . Distinguishing the both cases, we 
(D.9) obtain 

I 



2s 



b + c + d = A"'(0+)A'(0+) f f c- 2q " s f dt 2 e- q i t2 4s 2 - f dt 2 c- q 2 t2 (t 2 2 +4s 2 )- f dt 2 e-& 2 4sh 

J qi ,q 2 J s>0 J J J 

Lo 2s 



-A"'(0+)A'(0 + ) 



1 



'91,92 mmi + qi) 

This integral can be simplified through symmetrization. Using that 



/ 



+ 



we obtain 



The next diagram is e: 



qM(q 2 i + ql) qtqiiql + q 2 2 ) 



b + c + d = --A"'(0+)A'(0 + )/ 2 



4 4 

q\q\ 



r 2 



(D.ll) 



(D.12) 



(D.13) 



U-t 3 ) 



e = / f e-^ t ^«2 2 *2-9 3 2 t 3A " (uo _ u _ t2 _ ti))A ' (u _ t2 _, 
Jqi-m Jt 1 ,t 2 ,t 3 

= 1 f e-^-^-^s [A"'(0 + )A'(0+)(t 2 +t 1 )sgn(t 3 -t 2 ) + A"(0+) 2 (t 3 -t 2 )] . (D.14) 

J 91 ,92 J tl,*2,*3 

By symmetrizing in (2 <-» 3), the term proportional to t\ and the term proportional to (£3 — t 2 ) vanish. The remaining term can 
be written as 



e = ~ f J e-^ tl ^ t2 -^ t3 A"'(0+)A'(0 + )|t 2 -i 3 | 



(D.15) 



91,92 *l,t2,t3 



Making the same change of variables to u and s as for (D.ll), the integration over u, s and t\ can be performed in this order, 
distinguishing the cases u < and u > 0. Both cases give the same result for a total of 



1 



L,g 2 q\q\(q\ + 



e = -A"'(0+)A'(0+) / 
This contains the new integral (given regularized) 

1,1 ^ l qi , q2 (??+m 2 )(^+™ 2 ) 2 (q 2 + q 2 + 2 m 2) 



(D.16) 



(D.17) 
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It is related to h, see (A.26) and I a, see (A.18): 

I v +h=I A . (D.18) 
It is calculated in appendix E. The last diagram to be calculated is f : 

91,92 Jti,t-2,t 3 

= - f f B -fa-<fa-<&** [A"'(0+)A'(0 + )(t 2 + h) + A"(0+) 2 (ti + i 3 )] 

J gi, 92 Jti,t 2 ,t 3 

= -2A"'(0+)A'(0+)/a - 2A"(0+) 2 / A • (D.19) 
APPENDIX E: THE INTEGRAL /„ 

We have to calculate the integral defined as 

/?? ^ l qug2 (g 2 +m 2 )(g 2 + m 2 ) 2 ( g 2 + (Z 2 + 2m 2 ) • (E.1) 



This is done as follows: 



I = [ [ ^ e -a{ql+ m 2 )-[3{ql+7n 2 )- 1 (ql+ql+2ra 2 

J H ,92 J a,/3,7>0 

7-07 

\Jg / Ja,0, 

'ML 



2 

-m 2 (a+/3+2 7 ) 



Det a + 7 7 



-d/2 



/3e" 

/3, 7 >o \ 7 /? + 27^ 

/ 3 7 3-V ro ^ Q+ ' 3+2 )(l + 2a + /3 + a/3)- d / 2 

9 / «^a,/3,7>0 

, \ 2 

e" 9 ' T(4 - d) m- 2e J (E.2) 



with 



oo />oo 



J := / da/ d/3 — — — — — — — ¥ ? — —5—^ = Ji + J2 + J3 (E.3) 

J (1 + 2a + p + a/3) 2 (a + /? + 2) e 







/•OO />1 

Ji := / da/ d^y^-^ — f^- — ^ ^ ' o ^ ; ~ = 21n3 - 31n2 + 0(e) (E.4) 



/•oo /*oo 



J 2 := / da / d/3 



P (l + 2a + /3 + a/3 )2 

(1 + 2a + /3 + a/3) 2 (a + /3 + 2) £ 

/? (l + 2a + /? + a/?)4 



(1 + 2a + (3 + a/3) 2 (a + (3 + 2) £ (1 + a) 2 "f /3 1+ f 



= ln2-21n3 + 0(e) (E.5) 



1 2 

(1 + af-ip 1 



J 3 := I da I dp - ; 2 =- + 1 + Q(6) . (E.6) 



Thus 



' l! , 2 (?? + m2)(^ +m 2 )2(? 2 +g 2 +2m2) -( 2e2 + 4 € )( cJ i) 2 + Iinilc - (E7) 



91, 

r 



APPENDIX F: INTEGRALS IN LONG this is 

RANGE ELASTICITY CALCULATION 

i_ r ds ga r 

'(a) Jo s Jq 



(q 2 + m 2 )« r(c 

In the long-range case, there are two integrals which con- m - £ EiJ2 ^ [ c~ q2 ^\ (Fl) 

tribute to the renormalization of the disorder. At 1 -loop order T(a) \J ) 
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At 2-loop order, this is 



r(«) 



1 



1 ' !< (q 2 +m 2 )%(q%+m 2 ) a ((q 1 +q 2 ) 2 + m 2 )* 

91,92 

(F.2) 

This is evaluated as follows 

ds s§ f 00 df t Q f°° du it* 



r(«) 



a r(f ) 7 t i» 7 « r(f ) 

s(gj+m 2 )-i(^+m 2 )-ti(( 9 i+g 2 ) 2 +m 2 ) 



X / C 

•'91,92 



ds s 2 



oo df t a /-oo 



S r(f ) J t r(a) 7 « r(f ) 



du us 



det 



s + u u 
u t + u 



-{s+t+u)m 



(F.3) 



Making the replacement s — > su and t — > to and integrating 
over m, we obtain 



r(«) 



m 



r(c) 



r(a)r(f)2 



s,t>0 



( r+ fl +"t)° (rt+a+t) * (i+a+t) " 



r(e) 



r(a)r(f)2 



Jl 



(\ t r 

o Jo 



ds 



{Ji + J2 + J3 + 0(e)) 
(FA) 
(F.5) 



(st + s + t) a 

/OO /'OO 
dt / daa* -1 * 0-1 x 

[(si + s + i)- a -(l + s)- a t- a ] (F.6) 



/OO /'OO 
df / dsst-^i + s)!-"^- 1 - 

2 r(f)r(^) 



e r(a-f) 



(F.7) 



Ji and J 2 are now both integrated over s. Changing in J 2 the 
integration over t to that over 1/t, we obtain 



r(^) y 



df ■ 



(Ji + J 2 )| Q=1 = 

+ J 2)| Q=1 = 27rln2. 
Putting everything together, the final result is 



r(«) 



2e 2 + 4e Uo t(l + *)' 

, r'(a) r'(f) 



(F.8) 
(F.9) 
(F10) 



(F.ll) 



r(l) _ 
r(2) _ 



1 


+ 


In 2 


2T 2 " 


e 


1 


+ 


1 " 


2T 2 " 


4^ 



r(o) r(§) 
( e /^) 2 + o( e °) 



e/j Q) ) +0(e°) 
(F.12) 

(el^Y + 0(e a ) . (F.13) 



APPENDIX G: CALCULATION OF THE 
INTEGRAL /« 

The calculations for the corrections to friction are the same 
as with short-range elasticity, except that the integrals I\, I a 
and I v change. The first two have already been calculated in 

appendix F. We now attack the masterpiece, 1^ . For sim- 
plicity, we restrict ourselves to a = 1. 

7(D := / £ I i ^ 

' (q\+m 2 )2 (q 2 +m 2 ) [{q^+m 2 )^ + (g|+m 2 )2] 

(G.l) 



91,92 



Using: 



we have 



-= / dss- 3 / 2 c-^c- sx (G.2) 



1 = r 

y/a + Vb Jt 3 >o 

= _L y ( SlS2 )-ie"^ _ ^e- Slt 3 a - S2t 3 6 



t3,Sl,S 2 >0 



(G.3) 



With the help of (G.3), we can write 1^ as 



r(l) _ _1 L_ /" / +"^(' Sl s „\-| -*i(9?W)-(*2+si^)(92+™ 2 )-S2^(9 2 +m 2 ) -4^-4^ 

" ~47rr(i) J J 1 [12) 

91,92 tl,t2,t3,Sl,S 2 >0 



1 1 f f -q 2 Y f V(siS2)"5e V(fl+t2+(si+S2) '3>e"^^ 
_ 4^r(I)VV J 7 (iifc + sisrfi + (ai + s 2 )tit§ + s 2 *2*§) 1 -' 

tl,*2,*3,Sl,S2>0 
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1 1 



= ^W) il e ~ q ) 



-l+2e *1 2 (S1S 2 ) 2 e" 



*l,*2,*3 ? Si,S2>0 



(ii*2 + sis 2 + (si + s 2 )h + s 2 t 2 ) 1 2 

-1. , 3 ! 



ti,t 2 ,si,s 2 >0 



(*i*2 + sis 2 + (si + s 2 )ti + s 2 t 2 ) 1 2 (*i + t 2 + si + s 2 ) £ 



r(e) 



m J . 



(G.4) 



In the third line we have made the replacement ti — > and i 2 — > i^ 2 . In the fourth line we have integrated over i 3 . The 
integral J is again decomposed in converging parts (which can be evaluated at e = 0) and parts that can be integrated analytically: 



J = J1 + J2 + J3 + 0(e) 



Ji 
J2 
J 3 



h 2 (sis 2 ) 2 e 4s i 



*l,Sl,S2>0 



3--^- /* 

S1 4s 2 / 
Jti 



tl,Sl,S2>0 



/ 



*i 2 (sis 2 ) 2 e 4 =i~ 

_I _3 1_ 

h 2 («is 2 ) 2 c 4s i 



t 2 >l *1*2 + Sl«2 + Ol + S 2 )tl + s 2 t 2 t 2 (s 2 + ti) 

1 



o<t 2 <i *i*2 + sis 2 + (si + s 2 )ti + s 2 i 2 



t 2 3 (s 2 +ii) 



,si,s 2 >0 Jt 2 >l 

Integrating in J3 over i 2 , t\, si and s 2 (in this order) we find 



(G.5) 
(G.6) 

(G.7) 

(G.8) 

(G.9) 



In order to calculate J\ and J 2 , it is convenient to do the integration over t 2 in both integrals first. Taking the sum, some terms 
cancel: 



s_3 -J — 1 ln(s 2 +ii) -ln(sis 2 + ii(si +s 2 )) 
Ji + J 2 = / ^i 2 (sis 2 ) 2 c 4 -i 4 =2 — i '- i '-L 

/ti,si,s 2 >0 S2+I1 



(G.10) 



The logarithms have to be written as derivatives: 



T T 9 



f t^( SlS2 )-h ^-^((s 2 +t 1 )- 1+b -( S2 +t 1 )- 1 ( Sl s 2 +t 1 ( Sl + S2 )) b ) . (G.ll) 

',=0 Jt 1 ,S 1 ,S 2 >0 V ' 



Making the change of variables s\ — ► 1/si, s 2 — ► l/s 2 , ti — > ti/s 2 and si — ► sis 2 (in this order), we obtain 

d 



J1 + J2 



db 



[ *i h 4 b tf>((l + t 1 )- 1+b -(l + h)-hi b s 2 \l + h(l + s 1 )) b )e-i s ^ 1+s ^ . (G.12) 

=0 Jt 1 ,S 1 ,S 2 >0 V ' 



The integration over s 2 can now be done analytically: 



r(f-6) /i + Sl \ 6 -i r(|-26)(i + t 1 (i + Sl )) b /i + 



(l + ti) 1 -" 



6=0 ^tljSl 

In order to proceed, we split these integrals as follows 

Ji + J 2 = K 1 + K 2 + K 3 

d_ 

db 



(i + *i)*S 



Sl 



2b- 



K 2 



J 

-0 Jti,s 



s 1 i 1 



si 



Sl 2 i x 



(1 + ti) 1 - 6 
i r (f- 2& ) A 



96 



fc=0 Jti,si 



2b- 



K, = - — 



(i + ti)«! 
r(|)(i + i 1 (i + Sl )) b A + si 



fc=0 

d_ 

' db 



8^4- b r(i-&) 



Stt2 4- b r(l-26) 



(1+tl) 

r 



fc=0 

3 



(G.13) 

(G.14) 
(G.15) 

(G.16) 

(G.17) 



To evaluate if 3 one first has to take the derivative: 

ln(l+ti(l + «i)) 



K 3 



-8T 



(I) / 

Jti, 



Integrating first over t\ and then s\ gives 

2atanh(l/vTTsi) +ln(si) 



"1 v^rv^r(i+*i)(i + *i) ! 



K, 



= -47T 2 / 
J si 



(G18) 



sT (l + ai)* 
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= 8tt* (21n-27r) . 
Putting everything together his gives finally 



v 



1 

2? 



In 2 



(1)N2 



finite . 



(G.19) 



(G.20) 



APPENDIX H: FIXED-POINT FUNCTION AT 
SECOND ORDER 

In this appendix, we show how to obtain the fixed-point func- 
tion for A (it) at second order. We restrict the discussion to 
a = 2. We use the notations of equation (IV.5). First, one 
needs the 1-loop function y\(u) both by solving (IV. 11) nu- 
merically and as a Taylor-series about 0. The latter is obtained 
by deriving the 1-loop /3-function at the origin and fitting the 
coefficients as 



2/1 («) = 1_u +y 
139m 7 



u 
36 



u 

270 



u 



+ 



4320 17010 
571m 9 



204120 2351462400 
163879m 11 



5443200 

281m 10 
1515591000 

5221m 12 
354648294000 
5459m 14 

7447614174000 + 1830325659402240000 



2172751257600 

5246819m 13 
~ 10168475885568000 
534703531m 15 



(H.l) 



The /3-function at second order yields a linear differential 
equation for j/2 (u). It is numerically singular at small u. 
Therefore one has to expand it in a Taylor-series about 0. Us- 
ing the above information and the knowledge of £2, one finds 



both towards and towards infinity, one verifies that Taylor- 
expansion and numerically obtained curve coincide in their 
respective domain of validity. This is shown on figure H.l. 
One also verifies that the numerically obtained function con- 
verges to for large u, thus the exponent (2 obtained above is 
correct. 

It is a good question to ask for how large e the fixed-point 
function A(m) = §yi(it) + jgy2{u) might be a good approx- 
imation for the true disorder correlator. Let us note that if 
one demands that A(m) > 0, thus that forces be never anti- 
correlated, this is only satisfied if 



e < e r 



1.6 



(H.3) 




-0.5 



FIG. H.l: The fixed-point function of the RG-flow 3/2 (w) at sec- 
ond order in e. Upper curve: Numerical integration. Lower curve: 
Taylor-expansion about 0. 



APPENDIX I: STABILITY OF THE FIXED 
POINTS 



y 2 (u) = -1.14012 m + 0.967798 m 2 - 0.202495 m 3 



-0.019299 m 4 + 0.00259234 m 5 



0.0015302m° 



+0.000286423 m 7 - 6.25533 10" 6 m 8 



-0.0000206648 u 3 



6.48801 10~ 6 m 10 



-7.8566910" 



i 11 + 1.88404 IO^m 1 



+1.24668 10~ 7 m 13 



3.13093 10~ 8 m 14 



(H.2) 



We now consider the stability of the periodic fixed point given 
in (IV.33). Define K[f] as 

K[f] := Jim i [^(A*(«) + «/(«)) -/3(A*(u))] . (LI) 
The eigenfunctions and eigenvalues are 

K[f] = Xf . (1.2) 



The differential equation for 1/2 (m) is then solved numeri- 
cally, starting at u « 0.5. By integrating from that point 



We find the following solutions (with x = u(l — u) and nor- 
malized to /(0) = 1) 
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M — 


e i 






1 

— 1 






A 2 = 


— e — 


3 s ' 


h 


= i - 


(6 + Ae)x 




A 3 = 


-4e- 


-5e 2 , 


h 


= i - 


(15 + 20e)x + (45 + 85e)x 2 




A 4 = 


-25e 
3 


140e 2 
9 ' 


h 


= i - 


(28+ 23 £ £ )x+(% 6 + 2 ^)x 2 


,-4004 | 185402c \ 3 
V 9 81 ' 


A 5 = 


-14e 


- 35e 2 , 


h 


= i - 


(45 + 225e)a; + (585 + 4500e)a; 2 - 


(2925 + 110475e )x 3 + ("f + 4247556 )x 4 


A 6 = 


-21e 


- 66e 2 , 


h 


= i - 


(66 + 5l7e)x + (1320 + 16148e)a; 2 


(11220 + 169928e)a; 3 + (42636 + 3672 _ 944e )a; 4 



_ ( 298452 _|_ 28055588c ^5 



(1-3) 



This shows that apart from the constant mode (the shift) dis- 
cussed in the text, the fixed point is stable. 

APPENDIX J: CALCULATION OF 
CORRELATION FUNCTIONS 

In this Appendix we show how to compute a correlation func- 
tion in the renormalized theory. As an example we study the 
periodic case, i.e. we compute the amplitude Ad in (IV. 40). To 
do that we assume that we are exactly at the fixed point. 

The correlation function is time-independent, as was shown 
in Section (III D), and takes the scaling form: 

(u q u- q ) nl = ^-A*(0) m - d F d (^) , (J.l) 
el\ \7TlS 

where we have restored the factor previously absorbed in A. 
The scaling function is universal and satisfies F(0) = 1 since 
our calculation was performed at zero external momenta in 
presence of a mass and F(z) ~ B/z d at large z. In d = 4 
one has F^z) = 1/(1 + z 2 ) 2 . We want to obtain the scaling 
function to the next order; in particular to compute Ad we 
need B = 1 + be + 0(e 2 ). The universal amplitude reads: 



A d = 



2S d 



(2n) d el 



(J-2) 



= (l + 6e)(2 + C )(— + 10g 



0(e 3 ) 



which yields: 



1 21< a » + 3 + 66 2 
Ad = ~8 e+ 108 6 



(J.3) 



Computing b requires computing diagrams with external mo- 
mentum which we do now. Let us use straight perturbation 



theory with A , as in Section (IIIC 1). One has 

(q 2 + m 2 ) 2 (u q u_ q ) = A o (0) - A (0+) 2 /(<?) (J.4) 

(J-5) 



i(q) = I 
jp 



(p 2 + m 2 )((p + q) 2 + to 2 ) 



Let us reexpress this by the renormalized dimensionless dis- 
order given in (111.35) and (III. 18) at u = 0: 

A o (0) = to £ (A(0) + A'(0+) 2 to £ /(0)) . (J.6) 
This gives: 

(q 2 +m 2 ) 2 (u q u- q ) 

= to £ (A(0) - A'(0+) 2 to £ (/(<z) - 1(0))) 

= to £ -La*(0)(1 - e\m\I(q) - 1(0))) , (J.7) 
eli £+1 

where we have reestablished the factor A(u) = ^i-A*(u) 

and used the fixed point condition A*'(0+) 2 = eA*(0) This 
substitution acts as a counter-term which exactly subtract the 
divergence as it should. The result is finite. Using that: 

=ll c - ;s (P+9/2) 2 -t(p-'?/2) 2 -( S +t)m 2 
Jp Js,t>0 

= [ e- p2 [ (s + t)- d / 2 e- q2 ^- ( - s+t)m2 

Jp Js.t>0 

= 1°-' 

Jp 



m"T 2 



/ (1+*)- 
Jt>o 



d/2 



(!+<) + 



t q' 



-e/2 



1 + t TO 2 

One obtains the scaling function in the form (z — \q\/m): 



. (J.8) 



Fdiz) 



(1 + z 2 ) 2 
] 

+ 



At 



tz 2 



e/2 



(J-9) 



^ |l + | jT 1 ds In(l + z 2 (s S 2 ))| + 0(e 2 ) 
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1 



(1 + z 2 ) 2 

C^oo 1 



4 + z 2 



ln2-ln (2 + z 2 - + j + 0(e 



?0 



, L l + -[-2 + 21nz 



}+0(6 2 ) 



We want to match at large z: 

F d (z) = - i {l + be)z' = l(l + e (lnz + fe) + 0(e 2 )) (J.10) 

z z 

The above result yields 

b = -l. (J.ll) 

APPENDIX K: ANOMALOUS AND 
NON-ODD GRAPHS 

In this Appendix we write all anomalous 2-loop graphs con- 
tributing to the correction of a non-analytic disorder. In a first 
step we make no assumption and give their general expres- 
sions: Already at that stage some cancellations are apparent. 
In a second step we consider the limit v — > + at T = 0. 
We check all cancellations given in the text and show that no 
additional singularities occur. The multiplicity factors are in- 
cluded in the given expressions. Of course since we want only 
corrections to disorder we will give only the expressions when 
the separations of the times between the two external response 
fields are much larger than the separations within each con- 
nected component. If this were not the case, as is needed e.g. 
in the calculation of a 2-point correlation function to order 
A 3 , the above expressions should be reexamined separately. 
Equivalently, the expressions given here are correct only for 
u > and may become incorrect at u = 0. 

Graphs which are odd need not be considered (see main 
text). Each remaining graph, e.g. Cj is written in the shorthand 
notation form: 

graph Ci = [ [ F cCl . (K.1) 

Jx,y Jti>0 

The only anomalous non-vanishing graphs of class A are: 

Ryt\Ryt 2 Rx—yt 3 

(K.2) 

Cl = 2<A'K - u*_ t3 _ t2 _ t4 )AV_ t3 - «» t3 _ t2 _ tl )>A» 

(K.3) 

c 2 = 2(A( u * t4 - t£ ta _ t3 )A"« - u\_ tl ))A»(u) (K.4) 
c 4 = 2A>g - u*_ t2 _ t3 _ t4 )A'(ul t4 _ t2 ul t4 _ tl )A"(u) 

(K.5) 

c 5 = -2A(u* tl _ ta _ t3 - u*_ ti _ u )A"( U y - ul t2 _ ti )A"(u) 

(K.6) 

F e = Ff = R x - y t 1 R x -yt 2 Rxt 3 Ryt i (K.7) 

ei = A"(u){A'(u% 3 - u*_ tl _ u )A'(u\ - ul t3 _ t2 )) (K.8) 
f 2 = 2A"(u){A'(u% 3 - u*_ t3 _ t2 _ tl )A>(ul t2 _ t3 u*_J) . 

(K.9) 



For the graphs d one easily sees that the following relations 
are exact (with no other assumption than u 7^ 0): 

F c = Fa = R y t 1 Ryt 2 Rx-yt 3 Rxt i (K.10) 
d 2 +d 4 = (K.11) 
d 6 + d 8 =0. (K.12) 

The only anomalous non-vanishing graphs of class B are: 

Fk = F\ = R xtl Rxt 2 Ry~xt 3 Ry-xt l (K.13) 

kj = cA"( U )<A"( U - ta - u*_ tl )A(ul t2 _ t4 «» tl _ ts )) 

(K.14) 

h = -cA"(u){A"(u*_ t2 -u*_ ti )A( u y_ ti _ t4 -u»_ ti _ t3 )) 

(K.15) 

Fi = Rxt 1 Rxt 2 Ryt s Ryt 4 (K.16) 

i 3 = -A"(«)<A'K - ^ tl _ t2 )A'K - J) • 

(K.17) 

All graphs of class C exactly vanish. For instance: 

m 2 = A(0)A"> - U- t3 )A'(u) (K.18) 
n 2 = -A(0)A"> - «-t 3 )A'(u) . (K.19) 

We now evaluate these graphs in the quasi-static depinning 
limit, substituting A(u) by its power series as a function of u, 
as explained in the main text. We need in addition to (H.9): 

A» = A'(0+)sgn(u) + A"(0+)w + . . . 

A"(«) = 2A'(0 + )<5(u) + A"(0+) + . . . . (K.20) 

In A"(u) evaluated at zero we have written the <5-function 
which may in principle be needed. If this were the case that 
would pose two unpleasant problems: Firstly a different view- 
point were to argue that A" (u) should simply be continued to 
zero which does not pose any problem since it is pair. Sec- 
ond it would open the possibility to problematic singular terms 
(S(v) or 1/v) as v — > + . Fortunately, in all our 2-loop calcu- 
lations this never happens: these (5-functions, if put by hand, 
cancel. This confirms that, at least to this order, no pathology 
arises. 

Let us start with the sum c 2 + c 5 . Using (H.9) and (K.20) 
one sees that the term proportional to A(0)A"(0 + ) cancels. 
Let us test the <5-function. Then one needs to go one order fur- 
ther in the expansion of the A term since averages of the type 
5{u\)u2 have dimension one, similar to (sgn(wi)sgn(w 2 )), 
and can thus yield a non-zero result at zero temperature 
(higher order terms yielding dimensions as positive powers of 
the field are not needed as they vanish at zero T). This yields 
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c 2 + c 5 = 4A'(0+) 2 A»((K t4 - u*_ t2 _ t3 \ - \u*_ tl _ t2 _ t3 - ul ti _ ti \)6( U y - ul t2 _ tl )) , (K.21) 

which strictly vanish upon the replacement — u*, — ► v(t — t'). This is fortunate since this term would have led to a 1/v 
singularity. Note that all diagrams a — g in the 2-loop correction to r\ could a priori suffer from the same problem since A" 
functions must be expanded. However one notes that their arguments are always strictly positive in the depinning limit, which 
avoids, as it did here, the problem. Similarly one has 

c 4 = 2A'(0+) 2 A»(sgnK - u x _ t2 _ t3 _ u )sgn(ul t4 _ t2 u y _ u _ tl )) . (K.22) 



Performing the replacement — uf, — > v(t — t'), since the 
ti > and because F c is symmetric in t\ «-> ti one finds that 

c 2 + c 5 = c 4 = (K.23) 

at depinning. Note that these cancellations do not happen any 
longer, if the field is not a monotonic function, a question 
which will be discussed in Ref. 34 . 

A similar calculation shows that at depinning one has also: 

ki + li = . (K.24) 

There, in the singular part, the <5-function implies that t\ = t 2 



yielding the cancellation via a slightly different mechanism 
than above. 

Finally we are left with the only non-zero anomalous non- 
trivial graphs ci, e\, fi and i\ to compute, which is done in 
the text. 
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